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—'/Measurements  of  the  dependence  of  snow  albedo  on  wavelength,  zenith  angle,  grain  size,  impurity 
' content,  and  cloud  cover  can  be  interpreted  in  terms  of  single-scattering  and  multiple-scattering 
radiative  transfer  theory.  Ice  is  very  weakly  absorptive  in  the  visible  (minimum  absorption  at / = 0.46 
-pin)  but  has  strong  absorption  bands  in  the  near  infrared  (near  IR).  Snow  albedo  is  therefore  much 
lower  in  the  near  IR.  The  near-IR  solar  irradiance  thus  plays  an  important  role  in  snowmelt  and  in  the 
energy  balance  at  a snow  surface.  The  near-IR  albedo  is  very  sensitive  to  snow  grain  size  and 
moderately  sensitive  to  solar  zenith  angle.  The  visible  albedo  (for  pure  snow)  is  not  sensitive  to  these 
parameters  but  is  instead  affected  by  snowpack  thickness  and  parts-per-million  amounts  (or  less)  of 
impurities.  Grain  size  normally  increases  as  the  snow  ages,  causing  a reduction  in  albedo.  If  the  grain 
size  increases  as  a function  of  depth,  the  albedo  may  suffer  more  reduction  in  the  visible  or  in  the  near 
IR,  depending  on  the  rate  of  grain  size  increase.  The  presence  of  liquid  water  has  little  effect  per  se  on 
— snow^yrtical  properties  in  the  solar  spectrum,  in  contrast  to  its  enormous  effect  on  microwave 
emissiviTy?Snow  albedo  is  increased  at  all  wavelengths  as  the  solar  zenith  angle  increases  but  is  most 
sensitive  arountvX  = 1 -two-.  Many  apparently  conflicting  measurements  of  the  zenith  angle  dependence 
of  albedo  are  diffkiur  to  interpret  because  of  modeling  error,  instrument  error,  and  inadequate 
documentationoTgrain  size,  surface  roughness,  and  incident  radiation  spectrum.  Cloud  cover  affects 
snow  albedo  both  by  converting  direct  radiation  into  diffuse  radiation  and  also  by  altering  the  spectral 
distribution  of  the  radiation^luud  eovei  normally  causes  an  increase  in  spectrally  integrated  snow 
'mbedo.  Some  measurementsof  spectral  flux  extinction  in  snow  are  difficult  to  reconcile  with  the 
spectral  albedo  measurements.  The  bidirectional  reflectance  distribution  function  which  apportions 
the  reflected  solar  radiation  among  the  various  reflection  angler  must  be  known  in  order  to  interpret 
individual  sate.lite  measurements.  It  has  been  measured  at  the  snow  surface  and  at  the  top  of  the 
atmosphere,  but  its  dependence  on  wavelength,  snow  grain  size,  and  surface  roughness  is  still 
unknown.  Thermal  infrared  emissivity  of  snow  is  close  to  100%  but  is  a few  percent  lower  at  large 
viewing  angles  than  for  overhead  viewing.  It  is  very  insensitive  to  grain  size,  impurities,  snow  depth, 
liquid  water  content,  or  density.  Solar  reflectance  and  microwave  emissivity  are  both  sensitive  to 
various  of  these  snowpack  parameters.  However,  none  of  these  parameters  can  be  uniquely 
determined  by  satellite  measurements  at  a single  wavelength;  a multichannel  method  is  thus  necessary 
if  they  are  to  be  determined  by  remote  sensing. 
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An  understanding  of  the  reflection,  absorption,  and  trans- 
mission of  light  by  snow  is  important  for  two  general 
applications.  The  first  is  the  calculation  of  the  radiation 
budget  of  snowpacks  and  the  planetary  radiation  budget  over 
snow-covered  surfaces.  This  is  important  both  for  hydrolo- 
gy, because  radiation  is  usually  the  dominant  component  in 
the  surface  energy  budget  of  snow,  and  for  global  climate 
modeling.  The  second  application  .is  for  planning  the  remote 
sensing  of  snowpack  properties.  This  requires  modeling  of 
the  optical  properties  at  high  spectral  detail. 

Considerable  progress  has  recently  been  made  in  under- 
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standing  the  optical  properties  of  snow  in  the  solar  and 
infrared  regions  of  the  spectrum.  The  most  recent  review 
article  was  that  of  Mellor  [1977].  However,  the  modeling 
papers  of  Wiscombe  and  Warren  [1980a]  (hereafter  WWI) 
and  Warren  and  Wiscombe  [1980]  (hereafter  WWII)  included 
considerable  review  material.  They  reviewed  earlier  theoret- 
ical models,  snow  albedo  observations,  and  the  complex 
index  of  refraction  of  ice  and  reviewed  (where  necessary  to 
compare  with  model  results)  observations  of  the  dependence 
of  snow  albedo  on  wavelength,  grain  size  or  age,  liquid  water 
content,  solar  zenith  angle,  cloud  cover,  snowpack  thick- 
ness, and  snow  density. 

This  article  gives  a mere  thorough  review  of  observations 
and  modeling  than  did  WWI,  treats  flux  extinction  and 
bidirectional  reflectance  as  well  as  albedo,  but  especially 
reviews  the  considera  te  work  done  since  WWI  and  WWII 
were  written,  pointing  out  topics  for  further  research  where 
theories  and  observat  ons  are  lacking  or  in  conflict. 

This  review  of  optical  properties  is  limited  to  the  parts  of 
the  electromagnetic  spectrum  which  are  important  for  deter- 
mining the  climatic  role  of  snow  and  for  affecting  snowmelt. 
These  are  the  solar  (0.3  pm)  and  thermal  infrared  (5 

£.\  <40  nm)  wavelengths.  (Radiation  of  wavelength  shorter 
than  0.3  pm  is  absorbed  in  the  upper  atmosphere  and  does 
not  reach  the  surface.)  Other  parts  of  the  spectrum  (micro- 
waves)  will  be  mentioned  only  for  the  purpose  of  contrastive 
analysis  in  the  discussion  of  the  far-field  assumption  in 
scattering  theory  and  in  the  discussion  of  remote  sensing. 


B.  Definitions 

Definitions  for  reflectance  are  given  by  Siegel  and  Howell 
[1972,  pp.  47-88]  and  by  Nicodemus  et  al.  [1977]. 

The  reflected  radiation  is  not  perfectly  diffuse  but  is 
unevenly  distributed  among  the  reflection  angles  according 
to  th*  bidirectional  reflectance  distribution  function 
(BRDF).  This  function  R has  units  sr-1: 


R(6 b,  <t>\  X)  = 


dl(ff,  <*>',  X) 
Po dF(Oo,  <fo,  X) 


where  (flb,  <M  is  the  incident  (zenith,  azimuth)  angle,  po  = 
cos  6b.  (O'.  d>)  the  reflection  angle,  X the  wavelength,  F the 
incident  flux  (on  a surface  normal  to  the  beam),  and  / the 
reflected  intensity.  Unless  the  surface  has  azimuthally  de- 
pendent surface  features,  such  as  the  sas tragi  oriented  with 
their  long  axes  parallel  to  the  prevailing  wind  at  the  south 
pole  [ Carroll  and  Fitch,  1981],  the  dependence  of  R on  both 
and  4>  reduces  to  a dependence  only  on  the  relative 
azimuth  - d>- 

The  albedo  a,  is  the  ‘spectral  directional-hemispherical 
reflectance’;  it  is  the  integral  of  R over  all  reflection  angles: 

a,(6b,  X)  = |o'  ft’  dp.'  Jo  ’*(<*>.  O',  <t>\  X)  d# 

More  simply  stated,  the  albedo  is  just  the  upflux  divided  by 
the  downflux  at  a particular  wavelength,  usually  measured 
just  above  the  snow  surface. 

The  albedo  for  bemispherically  isotropic  incident  radiation 
is  the  diffuse  albedo  a/. 


a/X) 


-*£ 


Poo»(po.  X)  dfiQ 


In  general,  the  albedo  depends  on  the  distribution  ol  incident 
radiation  with  angle.  ^ 

The  spectrally  integrated  albedo  is  what  is  measi»re#by 
unfiltered  radiometers: 


d,l6 b)  = 


fa,(6b,  X)F  ] (0,  X)  dk 
SF  i(0,  X)dX 


where  F 1 (0,  X)  is  the  spectral  downflux  of  solar  radiation  at 
the  surface.  The  value  of  d thus  depends  not  only  on  the 
snow  properties  and  on  the  sun  angle  but  also  on  the 
atmospheric  composition  (water  vapor  content,  cloud  thick- 
ness, etc.),  which  affects  the  spectral  distribution  of  the 
sunlight.  ’ 

It  is  often  convenient  to  normalize  the  bidirectional  reflec- 
tance relative  to  the  albedo.  Thus  in  Figure  14  below  is 
plotted  not  R but  rather  the  anisotropic  reflectance  function 
/: 


M.  «\  *'“*»)  = 9,  v - <t>o)/as(8o)  (1) 


The  spectra]  emissivity  e(0,  X)  depends  on  emission  angle 
0 and  is  equal  to  the  absorptivity  or  the  coalbedo  [ 1 - as(6 b, 
X)]  by  KirchhofFs  law  [Siegel  and  Howcil,  1972]. 

Deep  in  a homogeneous  snowpack  (uniform  density  and 
grain  size  distribution  and  far  from  any  boundaries)  the 
spectral  flux  Fj(X)  is  attenuated  approximately  exponen- 
tially: 

Fl(X,z  + Az)  = Fi(X,z)e-^'A'- 
where  k,(X)  is  the  flux  extinction  coefficient: 


dlnFi(X)  -I  dr  j. (X) 



ks_i  is  often  reported  in  units  of  geometric  depth,  but  it  is 
better  expressed  in  units  of  liquid  equivalent  depth  in  order 
to  avoid  effects  of  snow  density  variation. 

The  important  quantities  for  calculating  snowmelt  are  the 
surface  albedo  and  emissivity.  The  important  quantities  for 
the  earth  radiation  budget  are  the  planetary  albedo  and  the  8- 
to  I2-/im  window  en  issivity.  None  of  these  are  measured  by 
narrow  field  of  view  satellites,  which  instead  measure  the 
planetary  bidirectional  reflectance  R over  a particular  wave- 
length band.  This  can  be  converted  to  planetary  albedo  for 
the  same  wavelength  band  if  the  anisotropic  reflectance 
function/is  known.  Further  conversion  to  a surface  albedo 
requires  knowledge  of  the  atmospheric  vertical  structure. 

Other  symbols  used  in  the  paper  arc  as  follows: 


g{x , m)  single-scattering  asymmetry  parameter; 
k,(X)  absorption  coefficient  of  pure,  bubble-free, 
polycrystalline  ice; 

k,(X)  flux  extinction  coefficient  for  snow; 
m(X)  = mn(\)  - imim(X)  complex  refractive  index  of  ice; 
g„,(x,  m)  single-scattering  extinction  efficiency; 
r snow  grain  radius; 
x size  parameter,  equal  to  2 w/X; 
ps  snow  density; 

<5  single-scattering  albedo. 

C.  Optical  Constants  of  Ice 

Theoretical  models  of  the  optical  properties  of  snow 
require  as  input  the  laboratory  measurements  of  the  refrac- 
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tive  index  mrt  and  absorption  coefficient  k,  of  pure  ice  as 
functions  of  wavelength.  They  are  combined  as  the  complex 
index  of  refraction  m = mn  - wn;m,  where  k,-  = AimtiJk. 

Of  particular  importance  for  solar  albedo  calculations  are 
the  recent  measurements  of  k,<A),  0.4  Aim  £ A £ 1.4  pm,  by 
Grenfell  and  Perovich  {1981],  The  absorption  coefficient  is 
so  small  in  the  visible  wavelengths  that,  to  measure  it 
accurately,  blocks  of  bubble-free  ice  as  long  as  2.8  m had  to 
be  grown  in  order  to  obtain  sufficient  light  attenuation. 

For  1.4  £ A £ 2.8  jtm  we  recommend  the  values  of  m(A) 
compiled  by  WWI,  and  for  2.8  £ A £ 33  #un  we  recommend 
the  measutements  by  Schaaf  and  Williams  [1973].  The 
optical  constants  of  ice  from  45-nm  to  8.6-m  wavelength  are 
reviewed  by  S.  G.  Warren  (unpublished  manuscript,  1981). 

D.  Measurements  of  Reflection  and  Transmission 
of  Light  by  Snow 

1 .  Albedo 

All-wave  albedo  has  been  routinely  measured  on  polar 
expeditions  for  many  years.  Time  series  of  albedo  show  high 
all-wave  albedos  (75-90%)  in  late  winter  and  early  spring, 
dropping  as  snowmelt  begins  to  about  60%.  Such  time  series 
have  been  reported  for  Greenland  by  Ambach  [1963]  and 
Diamond  and  G rdel  [1956];  for  Barrow  (Alaska)  by  Maykut 
and  Church  [1973];  for  McCall  Glacier  (Brooks  Range, 
Alaska)  by  Wendler  and  Weller  [1974];  for  McGill  Ice  Cap 
(Canada)  by  Havens  [1964];  and  for  snow-covered  sea  ice  in 
the  Antarctic  by  Weller  [1968]  and  in  the  Arctic  by  Langle- 
ben  [1971],  Summaries  of  monthly  average  or  seasonal 
average  all-wave  albedos  have  been  reported  for  Antarctic 
stations  by  W.  Schwerdtfeger  [1970,  p.  258]  and  for  drifting 
ice  islands  in  the  Arctic  by  Chernigovskii  [1963,  p.  269].  For 
the  Antarctic  Plateau,  where  snow  never  melts,  Schwerdt- 
feger found  that  ‘a  useful  value  of  0.8  as  the  lower  limit  of 
surface  albedo  appears  to  be  certain.’ 

Apart  from  these  routine  measurements,  problem-directed 
research  has  sought  to  identify  the  factors  influencing  snow 
albedo.  The  albedo  of  both  dry  snow  and  melting  snow  is 
normally  found  to  increase  as  solar  zenith  angle  increases,  as 
measured  by  Hubley  [1955],  Liljequist  [1956],  Rusin  [1961], 
Bryazgin  and  Koptev  [1969],  Korff  et  al.  [1974],  and  Carroll 
and  Filch  [1981].  These  measurements  are  examined  in 
section  J2  below.  Some  workers,  however,  found  the  oppo- 
site trend.  Havens  [1964]  reported  highest  albedos  at  mid- 
day, as  did  Kondratiev  et  al.  [1964]. 

Cloud  cover  affects  both  the  spectral  distribution  of  irradi- 
ance  and  the  effective  incident  zenith  angle.  It  normally 
causes  an  increase  in  all-wave  snow  albedo.  An  increase  of 
5-10%  relative  to  clear-sky  albedo  was  found  by  Liljequist 
[1956]  and  11%  by  Weller  [1968],  both  on  the  Antarctic 
coast,  and  5-7%  by  Hanson  [1960]  at  the  south  pole. 
However,  Carroll  and  Fitch  [1981]  have  now  found  cloud 
cover  to  reduce  albedo  at  the  south  pole;  this  can  be 
attributed  to  the  unusually  steep  dependence  of  albedo  on 
zenith  angle  which  they  find,  described  in  section  J2  below. 
Grenfell  et  al.  [1981]  found  snow  albedo  to  increase  with 
cloud  optical  thickness  at  a mid-latitude  site. 

The  effect  of  snow  thickness  on  the  albedo  of  a thin 
snowpack  over  a black  surface  was  investigated  by  Giddings 
and  LaChapelle  [1961]  for  monochromatic  light  at  A = 0.59 
pm.  They  found  the  snow  albedo  to  reach  within  3%  of  its 
asymptotic  value  at  a depth  of  R 'liquid  equivalent),  as 


did  O’Neill  and  Gray  [1973]  for  broadband-filtered  sunlight, 
0.3-1. 1 ftm.  However,  Warren  and  Wiscombe  [1980,  p. 
2742]  think  that  these  snow  samples  were  probably  some- 
what contaminated  and  that  pure  snow  would  require  4 times 
this  thickness  to  reach  within  3%  of  the  semi-infinite  albedo. 

The  reduction  of  albedo  due  to  snow  aging  has  been 
documented  for  visible  wavelengths  by  Holmgren  [1971]  and 
Grenfell  and  Maykut  [1977]  and  for  the  near  IR  by  O’Brien 
and  Munis  [1975].  Grenfell  et  al.  [1981]  studied  the  progress 
of  spectral  albedo  changes  (0.4  £ A £ 2.5  pm)  due  to  snow 
aging. 

Spectrally  detailed  measurements  are  necessary  for  an 
understanding  of  the  physical  processes  affecting  snow  albe- 
do. Many  of  these  measurements  were  reviewed  by  WWI 
and  WWII.  The  most  accurate  measurements  are  probably 
the  following.  Albedo  measurements  in  four  spectral  bands 
for  clean  Antarctic  snow  were  made  by  Liljequist  [1956]. 
High  spectral  resolution  albedo  measurements  were  report- 
ed for  0.4  £ A s 1.0  pm  on  the  Arctic  Ocean  by  Grenfell  and 
Maykut  [1977],  for  0.4  £ A £ 1 .5  pm  at  South  Pole  Station  by 
Kuhn  and  Siogas  [1978],  for  0.4  £ A £ 2.5  pm  in  the  Cascade 
Mountains  by  Grenfell  et  al.  [1981],  and  for  0.34  £ A £ 1.1 
pm  on  the  Great  Lakes  by  Bolsenga  [1981]. 

Grenfell  and  his  colleagues  have  developed  a portable  (16 
kg)  scanning  spectrophotometer  for  field  measurements.  The 
instrument  described  by  Roulet  et  al.  [1974]  was  useful  for 
0.4  £ A £ 1.0  pm;  it  has  been  improved  by  Grenfell  [1981] 
with  the  use  of  a circular  variable  nterference  filter  for  A < 
1.375  pun  (resolution  AA  £ 0.03  pm)  and  fixed  wavelength 
filters  to  extend  the  wavelength  range  out  to  2.45  A«n  (AA  •« 
0.1  pm). 

2.  Bidirectional  Reflectance 

The  bidirectional  reflectance  measurements  of  O’Brien 
and  Munis  [1975]  were  designed  principally  to  investigate 
the  spectral  dependence  of  reflectance  for  0.6  £ A £ 2.5  pm; 
only  a narrow  range  of  incidence  and  detector  angles  was 
employed.  Although  they  are  not  albedo  measurements, 
they  were  used  as  proxy  evidence  for  near-IR  spectral 
albedo  and  its  dependence  on  snow  age  by  the  albedo- 
modeling efforts  of  Choudhury  and  Chang  [1979a,  b]  and 
Wiscombe  and  Warren  [1980a]. 

Measurements  of  spectrally  integrated  bidirectional  reflec- 
tance over  a large  range  of  angles  were  made  by  Dirmhirn 
and  Eaton  [1975].  Measurements  over  a restricted  range  of 
angles  but  for  a variety  of  snow  types  were  reported  by 
Middleton  and  Mungall  [1952].  Section  L below  reviews 
these  as  well  as  the  aircraft  measurements  of  Griggs  and 
Marggraf  [1967]  and  Salomonson  and  Marlatt  [1968a,  b]. 
The  dependence  of  the  BRDF  on  wavelength,  grain  size,  and 
surface  irregularity  has  not  been  adequately  studied  either 
experimentally  or  theoretically. 

3.  Flux  Extinction 

The  monochromatic  flux  extinction  coefficient  k,(A)  de- 
creases rapidly  with  depth  near  the  surface  where  a signifi- 
cant fraction  of  the  upwelling  radiation  escapes  the  snow- 
pack.  Below  a few  centimeters  depth  the  effect  of  the  top 
boundary  is  no  longer  noticeable,  and  one  measures  an 
‘asymptotic’  monochromatic  flux  extinction  coefficient 
which  is  independent  of  depth  for  a homogeneous  snowpack. 
The  asymptotic  flux  extinction  coefficient  has  been  mea- 
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TABLE  1.  Models  for  the  Optical  Properties  of  Snow 


Grain 

Wavelength 

Anisotropic 

Input 

Size 

Dependence 

Scattering 

Sun  Angle 

Thin  Snow 

Reference 

Parameters 

Enters 

Examined 

Considered 

Dependence 

Treated 

Comments 

Dunkle  and  Bevons 

f,  m 

as  t 

yes 

yes 

for  diffuse  incidence  and  high 

11956] 

albedo 

Giddings  and  LaChapelle 

1,  mu,, 

as  / 

yes 

for  diffuse  incidence  and  high 

[1961] 

albedo 

Barkstrom  (1972) 

(O,0q 

yes 

must  be  tuned 

Barkstrom  and  Querfeld 

U>,  g,  00 

yes 

yes 

yes 

yes 

unrealistic  g 

[1975] 

Bohren  and  Barkstrom 

r,  m 

yes 

yes 

for  diffuse  incidence  and 

[1974] 

high  albedo 

Berger  [1979] 

m,  p, 

yes 

yes 

yes 

for  high  infrared  emissivity 

Choudhury  and  Chang 

r,  m,  p 

yes 

yes 

yes 

yes 

for  diffuse  incidence  and 

[1979 a,  6] 

albedo  »0.1 

Wiscombe  and  Warren 

r,  m,  6b 

yes 
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yes 
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used  in  this  paper 

[1980a] 

Choudhury  and  Chang 

r,  m,  6b.  s2 

yes 

yes 

yes 

yes 

‘surface  reflection’  included 

[1981] 

Symbols  used  are  as  follows:  g,  single-scattering  asymmetry  parameter;  /,  photon  mean  path  length  through  ice;  m = mle-  im,„,  complex 
index  of  refraction  of  ice;  r,  snow  grain  radius;  s2,  variance  of  surface  facet  slopes;  t,  ice  lamina  thickness;  p,  single-scattering  backscattered 
fraction;  6b,  solar  zenith  angle;  p„  snow  density;  and  d>,  single-scattering  albedo. 


sured  as  a function  of  wavelength  by  Liljequist  [1956]  and 
with  better  spectral  resolution  by  Grenfell  and  Maykut 
[1977]  and  Kuhn  and  Siogas  [1978],  These  relatively  mono- 
chromatic measurements  are  better  suited  to  testing  theoreti- 
cal models  than  are  all-wave  extinction  measurements  and 
are  discussed  below  in  section  12. 

Only  a few  such  monochromatic  measurements  of  flux 
extinction  have  been  reported.  Far  more  commonly  mea- 
sured is  the  attenuation  of  all-wave  solar  radiation  in  snow, 
which  has  been  reported  by  (among  others)  Ambach  and 
Habicht  [1962],  Ambach  [19631,  Weller  [1969],  and 
Schwerdtfeger  and  Weller  [1977].  A number  of  other  mea- 
surements were  reviewed  by  Mellor  [1977],  who  also  clearly 
explained  intuitively  the  fact  that  *,  decreases  as  grain  size 
increases.  Unlike  the  monochromatic  k,(X)  the  all-wave  k, 
does  not  quickly  reach  an  asymptote.  It  decreases  with 
depth  because  of  the  changing  spectral  composition  of 
sunlight  with  depth.  At  great  depth,  where  all  but  the  blue 
light  is  filtered  out  by  the  snow,  ks  will  reach  the  asymptotic 
value  corresponding  to  X = 0.46  pun. 

4.  Intensity  Extinction 

The  extinction  and  scattering  of  a directed  beam  of 
monochromatic  radiation  as  a function  of  angle  and  depth  in 
the  snowpack  has  been  studied  by  Ambach  and  his  co- 
workers  [e.g.,  Ott,  1974].  These  measurements  could  be 
useful  for  testing  future  models  which  may  attempt  to 
calculate  intensities  as  well  as  fluxes  within  the  snowpack. 

E.  Modeling  the  Optical  Properties  of  Snow 

Modeling  of  the  reflection  and  transmission  of  light  by 
snow  has  nearly  a 30-year  history.  A rather  oversimplify 
summary  is  given  in  Table  1 . The  early  models  of  Dunkle  and 
Bevans  [1956]  (hereafter  DB)  and  Giddings  and  LaChapelle 
[1961]  did  not  explicitly  compute  scattering  by  individual  ice 
grains  but  set  up  a two-stream  radiative  transfer  framework 
which  required  two  input  parameters.  These  two  parameters 
can  be  loosely  related  to  an  effective  grain  size  and  an 
absorption  coefficient,  but  they  are  normally  found  by  fitting 


experimental  data.  Because  these  models  are  computational- 
ly simple,  they  have  been  used  extensively  for  fitting  albedo 
and  flux  extinction  data  [e.g.,  Weller,  1969;  P.  Schwerdt- 
feger, 1969;  Bergen,  1970,  1971,  1975;  Schlatter,  1972; 
O'Neill  and  Gray,  1973].  However,  they  are  not  generally 
applicable  outside  'be  wavelength  range  where  they  are 
tuned. 

The  application  of  modem  radiative  transfer  theory  to 
snow  was  pioneered  by  Barkstrom  and  Bohren,  who  started 
with  tne  single  scattering  by  individual  ice  sphere  and  used  a 
number  of  approximations  to  relate  these  to  observable 
quantities.  Bohren  and  Barkstrom  [1974]  (hereafter  BB) 
obtained  very  simple  equations  which  (as  shown  below)  are 
applicable  only  for  the  visible  wavelengths.  Choudhury  and 
Chang  [1979a]  also  started  with  single  scattering  by  ice 
particles  and  used  a two-stream  method  for  radiative  trans- 
fer. They  did  not  attempt  to  derive  simple  parameterizations 
as  had  BB.  Their  model  is  applicable  over  a wider  wave- 
length range  than  is  BB’s  and  is  more  accurate  than  the 
model  of  DB. 

The  most  accurate  model  now  available  for  computing 
radiant  fluxes  in  snow  (short  of  a much  more  costly  doubling 
or  discrete  ordinates  method)  is  the  delta-Eddington/Mie 
theory  model  used  by  WWI.  Although  we  sometimes  use  it 
in  this  section  as  a benchmark  to  criticize  other  models,  one 
should  keep  in  mind  that  the  WWI  model  still  has  shortcom- 
ings, which  are  discussed  in  later  sections  of  this  paper:  it 
neglects  effects  due  to  close  packing  (which  restricts  its 
validity  to  X s 20  pm)  and  nonsphericity  of  snow  grains 
(which  may  cause  errors  at  very  large  solar  zenith  angles  fib 
-»  90°),  and  it  calculates  only  fluxes,  not  intensities,  so  it 
says  nothing  about  the  BRDF.  However,  it  has  proven  very 
useful  in  explaining  quantitatively  the  influence  of  snow 
parameters  and  environmental  parameters  on  spectral  albe- 
do. 

1.  Early  Two-Stream  Models 

Dunkle  and  Bevans  modeled  the  snowpack  as  a stack  of 
horizontal  ice  layers.  They  calculated  the  Fresnel  reflection 
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at  normal  incidence  on  each  layer,  as  well  as  the  absorption 
of  light  passing  vertically  through  the  layer  according  to 
Sauberer' s [1950]  measurements  of  Given  these  reflection 
and  absorption  coefficients,  as  a function  of  wavelength, 
they  used  the  Schuster  two-stream  method  to  examine  the 
dependence  of  albedo  on  the  thickness  of  the  ice  layers.  In 
Figure  1 we  reproduce  Figure  3 of  DB  together  with  the  more 
accurate  calculations  of  the  Wiscombe-Warren  model,  as- 
suming that  DB’s  layer  thickness  represents  the  snow  grain 
diameter.  The  DB  model  gives  tue  correct  qualitative  behav- 
ior of  afr,  X),  showing  the  decrease  of  albedo  due  to 
increased  grain  size  as  the  snow  ages,  as  well  as  the  fret  that 
albedo  is  iower  in  the  near  1R  than  in  the  visible.  However, 
there  are  substantial  errors.  These  errors  are  due  to  treating 
the  ice  as  sheets  rather  than  as  particles,  the  assumption  of 
normal  incidence  Fresnel  reflectivity,  and  the  use  of  two- 
stream  theory  (An  additional  small  part  of  the  discrepancy  is 
due  to  the  WWI  model’s  use  in  Figure  1 of  the  new 
measurements  of  m,m  by  Grenfell  and  Perovich  [1981].) 

DB  gave  formulas  for  both  albedo  and  transmittance  for 
both  thin  snow  and  semi-infinite  snow.  Although  they  ob- 
tained their  absorption  and  reflection  coefficients  frem  labo- 
ratory mesurements  on  pure  ice,  later  users  of  their  model 
have  treated  these  coefficients  as  two  adjustable  parameters 
to  be  fit  to  field  observations  of  snow. 

Giddings  and  LaChapelle  [1961]  (hereafter  GL)  used  a 
diffusion  model  which,  like  DB's,  also  employs  two  adjust- 
able parameters,  a diffusion  coefficient  and  an  absorption 
coefficient.  The  GL  model  is  actually  equivalent  to  *he  DB 
model,  because  the  diffusion  approximation  is  a form  of  two- 
stream  approximation.  (Recently,  both  Meador  and  Weaver 
[1980]  and  Zdunkowski  et  al.  [1980]  have  shown  that  all  two- 
stream  approximations  are  equivalent  and  can  be  put  into  a 
common  framework.)  The  diffusion  coefficient  was  related 
to  a length  I which  is  the  average  distance  a photon  travels 
through  ice  between  air-ice  interfaces,  so  it  is  interpreted  as 
the  effective  grain  diameter.  GL  estimated  that  the  simple 
diffusion  model  was  accurate  if  there  were  a large  number  of 
scatterings  before  the  absorption  of  a photon,  that  is,  if  ad  2- 

0.8.  Correcting  for  ’nomliffuseness’  (taking  into  account  the 
nonunit  ratio  of  downflux  to  upffux,  F]  IF  l < 1)  was 
thought  to  extend  the  validity  of  the  model  down  to  albedo  as 
low  as  0.5. 
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Fig.  1.  Model  calculations  of  semi-infinit;  diffuse  albedo  as  a 
function  of  wavelength  for  various  snow  grain  radii.  Dashed  lines 
are  calculations  by  Dunkle  and  Bevans  [1956,  Figure  3].  Solid  lines 
are  calculations  using  the  model  of  WWI,  with  the  new  nth,,  (X) 
measured  by  Grenfell  and  PerotUh  [1981], 


For  small  absorption  (large  albedo),  equations  13  and  21  of 
GL  imply  that  1 - ad  « rtn,  a result  also  obtained  later  by 
Bohren  and  Barkstrom  [1974]. 

GL  made  measurements  of  albedo  and  transmission  of 
thin  snow  over  a black  background  for  monochromatic  light 
(X  = 0.59  p m)  to  determine  the  two  free  parameters  in  the 
diffusion  model.  They  found  the  same  two  parameters  to  fit 
both  the  albedo  data  and  the  transmission  data.  However, 
the  value  of  I obtained  was  20  times  the  grain  radius 
estimated  by  eye,  I - 20r,  ^nereas  it  should  be  / = 2r 
according  to  the  interpretation  of  / as  grain  diameter. 

GL  also  analyzed  the  disti  -bance  of  the  radiation  field  due 
to  a radiometer  inserted  into  the  snow.  The  instrument 
measures  less  downflux  than  would  be  present  in  the  undis- 
turbed snow  because  it  is  blocking  some  of  the  upflux  that 
otherwise  could  be  scattered  back  down.  This  is  one  of  the 
reasons  why  transmission  measurements  are  more  difficult 
than  albedo  measurements. 

2.  Single-Scattering  by  Ice  Grains  Introduced 

The  scattering  and  absorption  of  radiation  by  a single  ice 
particle  are  described  by  three  quantities: 

1.  Extinction  efficiency  Qcxl  is  the  ratio  of  the  extinction 
cross  section  to  the  geometric  cross  section.  For  large 
particles  (r  » X),  Qex,  is  close  to  its  geometric  optics  limit  of 
2. 

2.  Single-scattering  albedo  to  is  the  ratio  of  scattering 
efficiency  to  extinction  efficiency.  It  is  the  probability  that  a 
photon  intercepted  by  a particle  will  be  scattered  rather  than 
absorbed. 

3.  Phase  function  P((l\,  ft2),  when  multiplied  by  to,  gives 
the  probability  that  a photon  incident  from  angle  ft,  = (0,, 
<Pi)  will  be  scattered  into  angle  ft2  = (02,  <fc)-  For  a spherical 
particle,  P is  a function  only  of  the  cosine  of  the  scattering 
angle  ft2  • ftt.  The  complete  phase  function  is  needed  for 
computing  intensity,  but  for  computing  fluxes  normally  only 
a single  measure  of  the  anisotropy  of  P is  needed,  commonly 
the  asymmetry  parameter  g,  which  is  the  mean  value  of 
fl2  - fi|,  or  the  backscattered  fraction  /3  [Wiscombe  and 
Grams,  1976;  Zdunkowski  et  al.,  1980]. 

Both  to  and  g are  dimensionless  with  ranges  0 £ to  ^ I and 
- 1 .s  g < 1 ; g = 0 corresponds  to  isotropic  scattering,  and  g 
= I to  completely  forward  directed  scattering. 

Barkstrom  [1972]  assumed  the  snowpack  to  be  semi- 
infinite, grey,  and  isotropically  scattering.  He  introduced  a 
zenith  angle  dependence  by  solving  the  radiative  transfer 
equation  for  intensity.  This  was  then  integrated  to  get  flux  in 
terms  of  the  X functions  of  radiative  transfer.  He  calculated 
that  albedo  would  increase  with  zenith  angle  in  approximate 
agreement  with  measurements  of  Rusin  [1961]  and  Liljequist 
[1956].  He  also  showed  that  the  (monochromatic)  flux  should 
decrease  faster  than  exponentially  at  the  surface  but  at  great 
depth  should  decrease  exponentially,  dF  i Idz  = -k^FI  , 
with  k,  independent  of  solar  zenith  angle. 

Hie  first  consideration  of  the  anisotropic  scattering  by  ice 
grains  was  that  of  Barkstrom  and  Querfeld  [1975],  who 
attempted  to  explain  the  bidirectional  reflectance  measure- 
ments of  snow  by  Middleton  and  Mungall  [1952].  Barkstrom 
and  Querfeld  used  the  adding-doubling  method  for  radiative 
transfer.  However,  in  order  to  match  Middleton  and  Mun- 
gall’s  measurements  they  required  quite  unrealistic  values  of 
asymmetry  parameter  (g  = 0.5,  corresponding  tor  = 0.1  pm. 
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whereas  g =>  0.9  for  a realistic  grain  size  r > 50  /xm). 

Bohren  and  Barkstrom  [1974]  used  geometrical  optics  in 
the  limit  of  small  absorption  to  calculate  the  scattering  by 
individual  ice  spheres.  They  obtained  g = 0.874,  close  to  the 
asymmetry  parameter  found  in  exact  Mie  calculations  for  a 
wide  range  of  particle  sizes  (g  = 0.89  for  X < 1 /xm  in  Figure 
4 of  WWI).  The  geometrical  optics  calculation  also  showed 
that  the  scattered  light  was  due  mostly  to  refraction  rather 
than  reflection. 

BB  made  a number  of  approximations  which  used  the 
assumption  that  c 1,  which  means  that  their  results 
apply  only  to  the  visible  wavelengths.  They  obtained  simple 
formulas  for  albedo  under  isotropic  illumination, 

ad  = 1 -8.43(Kfr),/2  (2) 

and  for  asymptotic  flux  extinction  coefficient, 

ks  = 0.65(k/t)i/2  (3) 

where  depth  is  measured  as  liquid-equivalent  depth.  Neither 
(2)  nor  (3)  involves  snow  density.  (BB’s  formula  did  show  ks 
to  be  proportional  to  snow  density,  but  this  dependence 
disappears  if  we  measure  depth  as  liquid-equivalent  depth. 
We  do  this  in  order  to  investigate  possible  near-field  effects 
which  are  ignored  in  all  published  models  and  which  would 
introduce  a density  dependence  into  (2)  and  (3),  as  discussed 
in  section  E5  below.) 

Figures  2 and  3 compare  equations  (2)  and  (3)  with  the 
more  accurate  results  of  WWI  for  several  wavelengths. 
Figure  2 shows  that  the  albedo  is  indeed  proportional  to  r'a 
for  X < 0.8  /xm.  Accordingly,  BB  found  good  agreement  of 
(2)  with  Liljequisl's  [1956]  observations  of  high  visible 
albedo  (ad  « 0.96)  using  Liljequist’s  measured  grain  size  r 
- 150  /xm;  this  is  the  wavelength  region  where  (2)  is 
applicable.  However,  (2)  becomes  useless  as  ice  becomes 
more  absorptive  in  the  near  IR.  At  X = 1 .3  /xm,  for  example, 
(2)  predicts  negative  albedos  for  r > 1 10  /xm.  Figure  3 shows 
that  the  flux  extinction  formula  (3)  is  better  behaved  than  the 
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Fig.  2.  Diffuse  albedo  versus  square  root  of  grain  radius  for  six 
discrete  wavelengths.  Solid  lines  are  calculated  using  WWI  model. 
Dashed  lines  are  calculated  using  equation  (42)  of  Bohren  and 
Barkstrom  [1974], 
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Fig.  3.  Asymptotic  flux  extinction  coefficient  k,  versus  (Mr)'12 
for  three  discrete  wavelengths.  Solid  lines  are  calculated  using  WWI 
model.  Dashed  lines  are  calculated  using  equation  (37)  of  Bohren 
and  Barkstrom  (1974). 


albedo  formula.  The  results  of  the  WWI  model  show  that  ks 
« r~'n,  as  (3)  predicts,  but  that  the  WWI  model  increasingly 
deviates  from  (3)  as  k,(X)  increases  (k,<1.5)  > ic,{1.8)  > k,<  1 .3 
/on)).  It  may  be  that  snow  albedo  and  flux  extinction  can  be 
parameterized  by  simple  formulas  like  (2)  and  (3),  but  more 
work  is  needed  to  develop  parameterizations  that  are  appli- 
cable over  wider  wavelength  ranges. 

Berger  [1979]  adapted  the  Bohren-Barkstrom  theory  for 
the  limit  of  large  absorption,  assuming  that  any  photon 
entering  an  ice  sphere  is  absorbed  by  it.  This  assumption 
makes  the  optical  properties  independent  of  grain  size. 
Berger’s  interest  was  to  model  the  infrared  emissivity  of 
snow,  and  his  large-absorption  approximation  is  reasonable 
for  r > 100  /um  in  the  thermal  infrared,  as  we  show  below  in 
section  M I • Berger  found  the  emissivity  e to  depend  on  snow 
density  p„  with  e increasing  as  p,  decreases,  owing  to  the 
reduced  average  angle  of  incidence  on  spheres  in  a regular 
array.  This  may  be  unrealistic,  because  the  derivation  de- 
pends on  the  particular  spherical  shape  of  ice  grains  and  on 
the  assumption  that  they  are  in  a regular  array.  However, 
Berger  found  the  dependence  of  e on  p,  to  be  weak  (Figure 
17  below). 

The  next  experimental  advances  which  stimulated  further 
modeling  were  the  spectral  bidirectional  reflectance  mea- 
surements of  O’Brien  and  Munis  [1975].  Choudhury  and 
Chang  [1979a,  b ] (hereafter  CCa  and  CCb)  used  the  Sagan- 
Pollack  two-stream  model,  which  was  rather  good  at  all 
wavelengths,  and  they  obtained  tolerable  agreement  with 
O’Brien  and  Munis’  measurements  (uncorrected  for  the 
reflectance  of  the  BaSO*  standard).  In  contrast  to  the  models 
of  DB,  GL,  and  BB,  none  of  which  were  applicable  for  ad  < 
0.5,  the  Choudhury-Chang  model  became  inaccurate  only 
for  flrfs0.I  (compare  Figure  4 of  CCa  with  Figure  1 of  CCb). 
In  their  two-stream  model,  CC  assumed  a backscatter  frac- 
tion (7.5%)  independent  of  wavelength,  using  single-scatter- 
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ing  albedo  calculated  from  r and  rn/m  by  means  of  a param- 
eterization due  to  Sagan  and  Pollack  [1967]:  <0  * } + j exp 
(- l)Ktr ),  where  b is  an  adjustable  parameter  taken  as  h ■ 1.67 
by  CCa  and  as  b = 2.0  by  CCb.  These  approximations  make 
the  CC  two-stream  model  more  accurate  than  the  DB  two- 
stream  model.  In  CCb  a ‘surface  reflection’  term  was 
introduced  which  calculated  the  Fresnel  reflection  from  a 
flat  sheet  of  ice  at  the  snow  surface.  However,  the  perfor- 
mances of  the  two  CC  models  were  not  compared  with  euch 
other  in  cither  of  these  papers,  so  the  effect  of  the  hypotheti- 
cal surface  layer  is  not  clear. 

3.  Wiscombe-  Warren  Model 

The  advances  in  modeling  made  by  WWI  were  to  use  Mie 
scattering  theory,  which  made  the  single-scattering  calcula- 
tions accurate  at  all  wavelengths  and  for  all  grain  sizes,  and 
to  use  the  delta-Eddington  approximation  [Joseph  et  al., 
1976]  to  handle  the  anisotropic  phase  function,  which  al- 
lowed the  model  to  calculate  albedo  for  any  solar  zenith 
angle  and  for  an  arbitrary  mix  of  diffuse  and  direct  radiation. 
The  detailed  measurements  of  visible  spectral  albedo  by 
Grenfell  and  Mayknt  [1977]  further  inspired  WWII  to  adapt 
their  model  to  calculate  the  effect  of  absorptive  impurities  on 
snow  albedo. 

The  snowpack  was  modeled  as  ice  spheres,  and  it  is 
argued  by  WWI  why  the  effects  of  nonsphericity  should  be 
small  in  relation  to  the  effects  of  grain  size  variation.  The 
scattering  and  absorption  of  light  by  single  ice  spheres  is 
described  by  Mie  theory.  Mie  calculations,  even  using  the 
fast  algorithms  of  Wiscombe  [1980],  are  extremely  time 
consuming  for  the  larger  snow  grain  sizes.  However,  for 
these  large  grains  (actually  for  any  grains  whose  size  param- 
eter is  x 2 100,  where  x is  the  ratio  of  the  circumference  of 
the  sphere  to  the  wavelength  of  light),  asymptotic  formulae 
have  been  developed  [Nttssenzvelg  and  Wiscombe,  19801 
which  arc  sufficiently  accurate  and  much  faster  than  the  Mie 
calculations. 

Mie  theory  assumes  that  the  particles  behave  as  isolated 
scatterers.  If  they  are  not  sufficiently  separated,  then  near- 
field effects  will  be  observed  that  arc  not  predicted  by  Mie 
theory.  WWI  examined  this  question  in  their  section  7 and 
concluded  that  the  near-field  effects  are  probably  negligible 
for  snow  in  the  solar  spectrum.  They  become  important  at 
longer  wavelengths,  and  a criterion  for  estimating  them  is 
described  below  in  section  E5. 

The  single-scattering  quantities  Qtx ,,  <i>,  and  g at  a particu- 
lar wavelength  are  functions  of  the  complex  refractive  index 
m and  the  effective  snow  grain  radius  r.  (This  is  the  area- 
weighted  mean  radius,  which  is  always  larger  than  the 
number-weighted  mean  radius.)  These  single-scattering 
quantities  become  the  input  to  a multiple-scattering  model. 
A logical  model  for  snow  albedo  is  the  delta-Eddington 
method,  beci'iu*.  it  can  adequately  handle  the  extreme 
asymmetry  of  scattering  by  ice  particles  in  snow,  in  which  a 
large  fraction  of  the  scattered  light  is  only  slightly  deflected. 
The  delta-Eddington  method  is  designed  to  be  efficient  and 
accurate  for  calculating  radiant  fluxes, 

The  model  can  also  be  used  to  calculate  snow  infrared 
emissivity  as  described  below  in  section  M. 

The  wavelength  dependence  of  snow  albedo  is  controlled 
by  the  variation  with  wavelength  of  the  absorption  coeffi- 
cient of  ice  k{\).  Within  that  constraint  the  model  shows  that 


snow  spectral  albedo  is  highly  sensitive  to  grain  size  and 
moderutely  sensitive  to  solar  zenith  angle  and,  in  the  visible 
wavelengths  only,  to  trace  amounts  of  absorptive  impurities. 
The  principal  results  of  the  model  are  summarized  in  the 
appropriate  sections  below. 

4.  Choudhury-Chang  Model 

Choudhiiry  and  Chqpg  [1981]  (hereafter  CC81)  and 
Choudhury  [1981]  have  now  abandoned  the  two-stream 
approach  in  favor  of  the  delta-Eddington  method.  ( Dozier  et 
al.  [1981]  have  shown  that  although  not  obvious,  equation 
(28)  of  CC81  is  indeed  equivalent  to  equation  (4)  of  WWI. 
This  is  the  special  case  of  diiect  incidence  on  a semi-infinite 
snowpack,  where  a*  is  a function  of  only  the  three  parame- 
ters w,  g,  and  Mo.) 

Instead  of  doing  Mie  calculations,  CC81  used  the  Sagan- 
Pollack  approximation  for  w mentioned  above  and  a param- 
eterization for  g which  they  devised  to  mimic  some  pub- 
lished results  from  Mie  theory.  However,  these 
approximations  are  quite  good  at  mimicking  the  Mie  results, 
at  least  for  the  larger  grain  sizes.  Thus  the  spectral  albedo 
calculations  of  CC8I  would  be  nearly  identical  to  those  of 
WWI,  except  for  the  use  by  CC8I  of  a special  ‘surface 
reflection*  term.  This  feature  of  the  CC8I  model  has  been 
criticized  by  Warren  and  Wiscombe  [1981],  who  think  that 
CC8!'s  special  accounting  of  surface  reflection  is  unneces- 
sary for  ordinary  snow  in  the  solar  spectrum.  Warren  and 
Wiscombe’s  main  points  are  as  follows:  (1)  The  nature  of  the 
interaction  of  electromagnetic  radiation  with  a snowpack 
depends  on  the  ratio  d/\,  where  d is  the  intcvparticlc  (center 
to  center)  separation;  snow  does  exhibit  a ‘surface*  to  radio 
waves  but  not  to  sunlight.  (2)  Even  if  one  does  want  to 
include  a surface  reflection  treatment  for  special  situations 
that  would  require  it,  it  is  formulated  incorrectly  by  CC8I. 
(3)  The  reason  for  introducing  surface  reflection,  namely,  to 
match  wavelength-integrated  albedo  observations,  is  insuffi- 
cient, because  the  model-measurement  discrepancy  could  in 
this  case  be  due  to  errors  in  the  atmospheric  radiation  mode! 
instead  of  the  snow  albedo  model. 

There  is  also  no  need  to  invoke  a special  surface  reflection 
to  explain  the  enhanced  specular  reflection  peak  at  low  sun 
angle.  The  explanation  for  specular  reflection  in  terms  of 
standard  single-scattering  and  multiple-scattering  theory  is 
included  below  in  the  section  on  zenith  angle  dependence 
(section  J). 

Although  CC8I’s  use  of  a surface  reflection  term  for  a 
homogeneous  snowpack  of  small,  randomly  oriented  grains 
seem.,  inappropriate,  such  a separate  modeling  of  surface 
reflection  could  indeed  be  tequired  for  a highly  nonrandom 
surface,  in  particular  for  the  case  of  glazed  crust,  or  ‘fim- 
spicgcl’  [LaChapelle,  1969,  Figures  59  and  60]. 

5.  Neglected  Effects 

a . Near-field  effects.  Because  snow  particles  are  close- 
ly packed,  they  may  be  in  each  other's  ‘near  field,'  meaning 
that  Mie  scattering  theory  is  inapplicable.  The  problem  of 
near-fleld  interference  was  mentioned  by  BB,  who  cited 
experiments  by  Blevln  and  Brawn  [1961]  on  the  density 
dependence  of  the  albedo  of  pigments  as  evidence  that  near- 
field effects  would  be  unimportant  for  • >w  of  p,  < 0.45  g 
cm-*.  But  in  these  experiments,  r **  X;  it  is  likely  that  near- 
fleld  effects  can  be  ignored  in  snow  up  to  considerably  higher 
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shall  thickness  or  wavelength  (for  snow  grain  radios  r * 200  flm) 


Fig.  4.  Transition  from  region  of  validity  of  far-field  approximation  to  near-fieid  regime.  The  assumption  is  made 
that  single-scattering  properties  of  a sphere  are  influenced  only  by  that  part  of  the  surrounding  medium  which  is  within 
one  wavelength  distance  of  the  surface  of  the  sphere.  A spherical  shell  of  thickness  X.  (surrounding  the  ice  particle  of 
radius  r)  contains  mostly  the  matrix  material  (air)  if  the  shell  is  thin  (X  < r)  but  contains  an  increasing  fraction  of  ice  from 
other  spheres  as  the  shell  expands.  The  ratio  of  the  density  in  the  shell  to  the  density  of  the  bulk  medium  (spheres  plus 
air)  is  plotted  as  a function  of  size  parameter  or  shell  thickness  for  five  different  regular  arrangements  of  spheres:  (1) 
hexagonal  close  packing  (ftooW  = 0.68  g/cm3),  (2)  simple  cubic  packing,  spheres  in  contact  'p5novi.  = 0.48),  (3)  simple 
cubic  packing,  spheres  not  quite  in  contact,  interparticle  gap  of  0.22/  (p^.  = 0.35)  (4)  interparticle  gap  of  0.56r  (psno»  = 
0.23),  (5)  interparticle  gap  of  1.33/  = 0.10). 


density.  WWI  reviewed  the  possible  near-field  effects  and 
pointed  out  that  interparticle  interference  should  be  neglect- 
ed for  particles  whose  center-to-center  separation  d is  large 
in  comparison  to  the  wavelength  X.  Since  d s>  X it:  the  solar 
spectrum,  no  interference  should  be  observed,  and  this  is 
confirmed  by  Bohren  and  Beschta' s [1979]  observation  that 
the  albedo  of  a thick  snowpack  is  independent  of  density. 
For  microwaves,  where  d s X,  interference  effects  will  arise, 
making  snow  microwave  emissivity  a function  of  snow 
density.  Both  <5  and  g are  altered.  It  is  possible  to  make  a 
rough  estimate  of  the  effect  using  Gate's  [1973]  adjustment 
to  Mie  theory  to  investigate  this  dependence  on  density. 
Gate  simply  altered  mrr  of  the  medium  to  be  not  that  of  air 
but  rather  a volume-weighted  mean  of  m„  (air)  and  m„  (ice). 
The  ‘medium’  should  probably  be  taken  to  be  a shell,  one  or 
a few  wavelengths  thick,  surrounding  the  particle.  The 
transition  from  the  far-field  regime  to  the  near-field  regime 
should  then  be  described  qualitatively  by  Figure  4.  Roughly 
stated,  snow  is  safely  in  the  far-field  regime  for  X < 1 pm  and 
in  the  near-field  regime  for  X > 1 cm,  with  a transition  region 
whose  location  depends  on  snow  density. 

It  is  possible  that  near-field  effects  could  become  impor- 
tant for  flux  extinction  (and  for  the  albedo  of  thin  snow)  even 
when  they  do  not  affect  the  albedo  of  deep  snow.  This  is 
because  the  albedo  of  deep  snow  depends  on  g and  <5  but  not 
on  2«:>  whereas  the  asymptotic  flux  extinction  coefficient 
depends  on  £>„t  as  well.  Whether  the  near-field  approach  of 
snow  grains  could  in  some  wavelength  region  significantly 
affect  Qt jtt  but  not  g or  <5  is  an  open  question. 

b.  Surface  irregularity.  The  model  calculations  of  albe- 
do all  assume  that  the  snow  surface  is  flat  Surface  irregular- 
ities can  reduce  the  albedo  relative  to  that  of  a flat  surface  if 


the  height  scale  of  the  irregularit.es  is  comparable  to  or 
larger  than  the  length  scale  and  if  both  scales  are  comparable 
to  or  larger  than  the  penetration  depth  of  light  (which 
depends  on  p„).  This  means  that  surface  roughness  of  210- 
cm  amplitude  (such  as  suncups)  is  required  to  reduce  visible 
albedo  but  that  much  smaller  irregularities  can  affect  near-IR 
albedo.  The  albedo  is  reduced  relative  to  that  of  a flat  surface 
because  some  of  the  light  reflected  from  a slope  does  not 
escape  to  space  but  instead  is  intercepted  by  the  slope  facing 
it.  There  is,  of  course,  no  enhancement  of  absorption  if  the 
albedo  of  the  flat  surface  is  0.0  or  1 .0;  the  enhancement  must 
reach  a maximum  at  some  intermediate  value  of  albedo. 
Pfeffer  [1982]  has  developed  a model  for  use  in  calculating 
the  enhanced  absorption  due  to  glacier  crevassing,  and  his 
model  might  be  adapted  to  the  study  of  smaller  irregularities. 

F.  Effect  of  Snow  Grain  Size  on  Albedo 

Figure  1 shows  the  calculated  spectral  albedo  of  snow  for 
diffuse  incident  radiation  for  three  different  grain  sizes 
expressed  as  spherical  radii  r.  These  radii  were  chosen  for 
comparison  with  Figure  3 of  DB,  but  the  smallest  size  is 
probably  unrealistic.  To  match  reflectance  measurements  of 
O’Brien  and  Munis  [1975]  for  new  snow,  WWI  never  needed 
grain  sizes  smaller  than  r = 50  pm.  What  the  optically 
equivalent  sphere  would  be  for  a nonspherical  snow  particle 
is  discussed  below. 

The  albedo  is  very  high  in  the  visible  wavelengths,  corre- 
sponding to  the  minimum  in  and  lower  in  the  near 
infrared.  The  albedo  drops  at  all  wavelengths  as  the  grain 
size  increases.  It  is  easy  to  understand  why  albedo  should 
decrease  with  increasing  grain  size.  Roughly  stated,  a pho- 
ton has  a chance  to  be  scattered  (or  a ray  to  be  bent)  when  it 
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crosses  an  air-ico  interface.  It  has  a chance  of  being  ab- 
sorbed only  while  it  is  passing  through  the  ice.  An  increase  in 
grain  size  causes  an  increase  in  the  path  length  that  musi  be 
traveled  through  the  ice  between  scattering  opportunities.  (A 
similar  dependence  of  reflectance  on  sand  grain  size  was 
lound  by  Leu  [1977]  for  beach  sand,  which  is  in  some  ways 
analogous  to  snow.) 

On  the  basis  of  matching  the  model  results  to  observations 
of  O'Brien  anil  Munis  [1975]  (see  Figures  10  orcl  I?  of  WWI) 
the  optical  grain  size  of  the  snow  surface  generally  varies  m 
the  range  from  50  mhi  for  new  snow  to  1 mtn  for  old  melting 
snow.  This  increase  of  grain  size  with  increasing  age  is  the 
normal  situation,  but  there  can  be  curious  exceptions.  Lilje- 
quist  [1956]  routinely  observed  the  albedo  at  Maudheim  (on 
a small  ice  shelf  at  the  coast  of  Antarctica)  to  he  low  after  a 
new  snowfall  and  to  rise  after  a windstorm.  At  the  very  coal 
temperatures  of  Maudheim,  snow  metamorphism  apparently 
proceeds  so  slowly  that  grain  size  changes  may  be  due  to 
other  effects.  Wind  probably  caused  a reduction  in  grain  size 
by  breaking  the  crystals  and  possibly  also  by  gravitational 
sorting.  The  latter  sould  cause  the  smallest  particles  to  settle 
out  last,  so  they  would  end  up  at  the  surface,  where  they 
would  dominate  the  albedo. 

1 . The  Optically  Equivalent  Sphere 

How  a field  measurement  of  snow  grain  size  translates 
into  the  radius  of  the  optically  equivalent  sphere  is  a subject 
of  current  research.  O'Brien  and  Koh  [1981]  found  that  the 
‘equal  projected  area’  assumption  gave  an  overestimate  of 
the  radius  of  the  equivalent  sphere.  Comparison  of  Grenfell 
et  a/.’s  [1981]  albedo  measurements  with  Figure  1 shows  that 
their  report  of  grain  radius  as  half  the  minimum  dimension  of 
the  snow  grain  is  a factor  of  2-5  smaller  than  the  optical  grain 
size.  The  best  conversion  procedure  must  lie  between  these 
two  extremes  and  is  probably  that  suggested  by  Dobbins  and 
Jizmagian  [1966]  (and  by  WWI).  Namely,  the  optically 
equivalent  sphere  is  that  which  has  the  same  volume-to- 
surface  ratio  V/S  as  the  nonspherical  snow  particle.  The 
results  of  Pollack  and  Cuzzi  [1980]  also  suggest  this.  (In  the 
case  of  a sphere,  of  course,  all  three  definitions  of  grain  size 
converge.) 

However,  any  attempt  to  treat  a nonspherical  particle  as 
an  equivalent  sphere  involves  a compromise,  because  the 
sphere  with  the  correct  w may  not  have  the  correct  g. 
Furthermore,  the  sphere  which  gives  the  correct  w at  one 
wavele  gth  may  not  do  so  at  another  wavelength.  The 
sphere  with  the  same  V/S  is  likely  to  be  most  appropriate  at 
wavelengths  where  absorption  is  small,  ntr  « I. 

2.  Grain  Size  Increasing  with  Depth 

The  albedos  in  Figure  1 are  for  a homogeneous  snowpack, 
that  is,  one  whose  average  grain  size  does  not  vary  with 
depth.  Grain  size  is  observed  to  increase  with  depth  in  a 
predictable  manner  in  the  Antarctic  but  can,  of  course,  both 
increase  and  decrease  with  depth  at  mid-altitudes.  In  order 
to  affect  the  shape  of  the  spectral  albedo  curve  the  grain  size 
must  change  rapidly  with  depth;  the  most  likely  situation 
where  this  effect  would  be  noticed  is  for  a very  thin  layer  of 
new  snow  on  a thick  layer  of  old  snow.  This  effect  has  been 
calculated  by  S.  0.  Warren  and  W.  J.  Wiscombe  (unpub- 
lished data,  1981).  An  increase  of  r with  depth  corresponding 
to  about  one-third  the  rate  of  increase  found  in  the  top 
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centimeter  at  the  south  pole  [ Stephenson . 1967,  Figure  8], 
causes  a greater  decrease  in  visible  albedo  than  in  near-IR 
albedo.  This  is  due  to  the  greater  penetration  depth  of  the 
visible  light,  which  thus  ‘sees’  a larger  average  grain  size. 
For  faster  rates  of  radius  increase  with  depth,  however, 
another  effect  takes  over.  The  albedo  drops  more  in  the  I- 
fj.m  region  than  in  the  0.4-^an  region.  Reference  to  Figure  I 
explains  this:  the  albedo  is  insensitive  to  grain  size  where 
albedo  is  very  high  (0.4  pm)  or  very  low  (2.8  *tm)  but  highly 
sensitive  in  the  region  of  intermediate  albedo  (I 

3.  Density 

There  have  been  many  reports  of  snow  albedo  decreasing 
as  density  increases.  Neither  BB  nor  WWI  nor  CC  obtained 
a density  dependence  in  their  models.  The  observed  depen- 
dence of  albedo  on  density  might  actually  be  a dependence 
on  grain  size,  since  density  normally  increases  as  grain  size 
increases.  Bohren  and  Beschta  [1979]  isolated  the  two 
parameters,  finding  albedo  unchanged  when  density  was 
artificially  increased  at  presumably  constant  grain  size. 

Density  enters  Bergen’s  [1975]  model  as  a parameter  used 
in  computing  the  air  permeability,  which  is  used  to  compute 
V/.5,  which  in  turn  is  used  to  compute  the  reflection  coeffi- 
cient in  the  DB  model.  We  interpreted  grain  size  above  as 
proportional  to  VIS , so  it  appears  that  Bergen’s  dependence 
of  albedo  on  density  could  actually  be  translated  into  a grain 
size  dependence,  when  grain  size  is  defined  as  above. 

G.  Effect  of  Liquid  Water  Content 

WWI  cited  both  experimental  and  theoretical  evidence 
that  the  effect  of  liquid  water  on  snow  albedo  is  simply  to 
increase  the  effective  grain  size,  because  the  refractive  index 
contrast  between  water  and  icc  is  very  small.  We  wish  only 
to  add  a footnote  to  that  statement  here.  O'Brien  and  Koh 
[1981]  have  pointed  out  that  the  slight  differences  (a  few 
percent  at  some  near-infrared  wavelengths)  in  reflectance 
noted  between  wet  melting  snow  and  the  subsequently 
refrozen  snow  (Figure  5)  are  in  the  right  direction  and  of  the 
right  magnitude  to  be  attributed  to  the  difference  in  spectral 
mim  of  water  and  ice.  This  is  most  obvious  between  1 .2-  and 
1 .4-/um  wavelength. 

H.  Effect  of  impurities  on  Snow  Albedo 

WWII  modeled  an  impure  snowpack  as  a mixture  of  ice 
particles  and  dust  or  soot  particles.  They  showed  that  small 
amounts  of  impurities  affect  snow  albedo  only  in  the  spectral 
region  where  absorption  of  light  by  ice  is  weakest,  mainly  in 
the  visible  (X  < 1 pm).  Reductions  of  visible  albedo  by  a few 
percent  can  be  caused  by  ~10  parts  per  million  by  weight 
(ppm  w)  of  desert  dust  or  ~0. 1 ppmw  of  carbon  soot.  A given 
amount  of  an  absorptive  impurity  causes  a greater  reduction 
in  albedo  for  coarse-grained  snow  than  for  fine-grained 
snow,  as  will  be  illustrated  in  Figure  7. 

1.  Soot 

WWI  and  WWII  showed  that  some  published  low  values 
of  visible  spectra]  albedo  for  apparently  dean  snow  could  be 
explained  by  the  model  only  if  the  measured  snow  sample 
had  contained  a grey  absorber  such  af  soot.  Discrepancies 
between  model  and  observation  were  obvious  for  the  recent 
careful  spectral  albedo  measurements  of  GwtfWf  .Way- 
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Fig.  5.  Effect  of  liquid  water  on  snow  spectral  reflectance.  (Figure  6 of  O'Brien  and  Munis  [1975].) 


kut  [1977]  at  Arctic  ice  island  T-3  and  Kuhn  and  Siogas 
[1978]  at  the  south  pole. 

Figure  6 shows  that  the  model  for  pure  snow  may  predict 
albedos  up  to  15%  higher  than  those  which  are  actually 
observed.  Unrealistically  large  grain  sizes  could  reduce  the 
calculated  visible  albedo  but  would  destroy  the  agreement 
with  observation  in  the  near  IR.  The  two  major  effects  which 
are  neglected  in  the  model,  namely,  nonsphericity  of  snow 
grams  and  near-field  scattering,  were  judged  (section  7 of 
WWI)  to  be  orders  of  magnitude  too  small  to  be  responsible 
for  the  discrepancy.  Insufficient  snow  depth  and  error  in 
laboratory-measured  ice  absorption  coefficient  were  also 
ruled  out  as  the  explanation,  leaving  impurities  as  the 
remaining  hypothesis.  In  order  to  match  the  spectral  shape 
of  albedo  at  T-3  (Figure  6)  a grey  absorber  such  as  soot  was 
implicated,  and  desert  dust  was  ruled  out. 

The  soct  in  snow  at  T-3  may  result  from  pollution,  either 
local  or  distant.  There  is  soot  in  the  Arctic  air  [Rosen  et  al., 
1981]  which  may  come  from  industrial  souices  in  Europe 
[Rahn,  1981].  but  it  is  questionable  whether  the  high  soot 
concentrations  found  necessary  to  explain  snow  albedo  at  T- 
3 (--0.2  ppmw  soot)  are  representative  of  the  entire  Arctic. 
H.  Rosen  (personal  communication,  1981)  found  soot 
amounts  in  the  range  0.01-0.06  ppmw  in  a preliminary 
experiment  on  snow  samples  from  Barrow,  Alaska. 

We  can  be  fairly  certain  that  any  soot  at  the  south  pole 
would  be  the  result  of  local  contamination  from  the  perma- 
nent camp.  For  representative  visible  spectral  albedos  for 
the  Antarctic,  one  should  therefore  use  Liljequist' s [1956, 
Figure  45]  high  values,  which  agree  with  the  pure-snow 
model,  rather  than  Kuhn  and  Siogas'  low  values.  For  X > 1 
pm,  soot  has  no  effect  (Figure  1 of  WWII),  and  Kuhn  and 
Siogas’  measurements  for  these  wavelengths  should  be 
representative  of  Antarctic  snow.  They  correspond  to  an 


average  grain  radius  of  100  pm,  which  agrees  fairly  well  with 
Stephenson's  [1967,  Figure  8]  grain  size  measurements  at 
Southice. 

Grenfell  et  al.  [1981]  have  recently  made  simultaneous 
measurements  of  spectral  albedo,  grain  size,  and  soot  for  a 
snowpack  in  the  Cascade  Mountains.  We  can  compare  these 
results  with  the  model  of  WWII.  The  measured  grain  size  in 
this  case  is  much  smaller  than  the  effective  optical  grain  size, 
as  mentioned  above,  so  we  ignore  the  grain  size  data  and 
instead  determine  the  optica]  grain  size  from  the  albedo 
measurements  at  X > 1 #un,  where  soot  has  no  effect  on 
albedo.  It  then  appears  that  in  order  to  explain  the  visible 
albedos  the  WWII  model  would  need  2-5  times  as  much  soot 
as  was  actually  found  in  the  snow.  A factor  of  2 difference 
can  be  explained  by  reconciling  the  definitions  of  soot 
amount.  Crhe  effect  of  a given  weight  fraction  of  soot  on 
snow  albedo  depends  on  the  assumed  values  of  absorption 
coefficient  and  density  for  soot.  Grenfell  et  al.  used  an 
operational  definition  of  soot  concentration  which  assumes  a 
mass  absorption  coefficient  of  8 m2  g_1  at  X = 0.55  pm.  In 
fact,  all  of  the  graphs  in  WWII  for  soct-containing  snow  are 
labeled  with  soot  amounts  that  are  probably  a factor  of  2 too 
large.  This  was  due  to  an  overestimate  of  soot  density,  which 
was  pointed  out  in  footnote  3 of  WWII  (see  also  Roessler 
and  Faxeot,  [1979]).  Figure  6 of  this  review  is  taken  from 
WW?.I  but  with  the  soot  amounts  now  altered  to  the  more 
likely  values.)  The  remaining  discrepancy  may  be  due  to  five 
causes. 

1 . The  new  accurate  laboratory  measurements  of  mim  by 
Grenfell  and  Perovich  [1981]  differ  somewhat  from  the 
values  of  Saaberer  [1950]  used  by  WWI  and  WWII.  Their 
use  causes  no  difference  in  the  maximum  albedo  value 
(although  its  position  is  at  0.46  pm  rather  than  0.40  pm)  but 
causes  a rise  of  ~-0.01  in  albedo  at  0.9  pm. 
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2.  Agreement  between  measurement  and  calculation  can 
be  further  improved  by  speculating  that  grain  size  increases 
with  depth  (which  was  observed  in  one  case). 

3.  Soot  which  is  reported  as  an  average  concentration  in 
the  snow  may  actually  be  concentrated  at  the  surface  (which 
was  also  observed  in  one  case)  where  it  has  more  effect  on 
albedo. 

4.  The  measurement  of  soot  concentration  may  have 
been  in  error.  Grenfell  et  al.  estimated  an  uncertainty  of  a 
factor  of  2-3  in  soot  concentration. 

5.  There  may  be  an  incorrect  assumption  in  the  model- 
ing. The  calculations  assumed  that  both  soot  particles  and 
ice  particles  were  surrounded  by  air.  The  possible  location  of 
soot  particles  inside  the  ice  grains  might  enhance  their 
absorption  of  light.  However,  it  has  not  been  demonstrated 
that  a significant  fraction  of  the  soot  could  be  inside  the  ice 
grains.  Soot  could  be  located  inside  a snow  particle  if  it  were 
attached  to  a dust  particle  which  served  as  an  ice  nucleus. 
But  it  seems  unlikely  that  soot  collected  by  scavenging  or 
dry  fallout  would  end  up  inside  ice  grains. 

2.  Volcanic  Ash 

Volcanic  ash  currently  seems  to  affect  snow  albedo  only 
locally,  but  at  times  in  the  past  it  has  probably  reduced  the 
albedo  of  the  entire  Antarctic  continent  (Gow  and  William- 
son [1971J,  reviewed  by  WWII). 

The  most  common  magma  is  basalt,  and  the  most  common 
volcanic  ash  is  andesite  (R.  Cadle,  personal  communication, 
1980).  Both  of  these  rocks  have  very  similar  optical  proper- 
ties for  short  waves:  m,„  ~ 1 x 10~\  constant  across  the 
visible  spectrum  [Pollack  et  al.,  1973].  Thus  their  effect  on 
snow  albedo  is  qualitatively  similar  to  that  of  soot,  that  is, 
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Fig.  6.  Comparison  of  calculated  snow  albedo  with  observations 
at  visible  wavelengths.  Solid  lines  are  measutements  of  Grenfell  and 
Maykut  (1977,  Figure  l]  made  in  summer  1974  at  ice  island  T-3  in  the 
Arctic  Ocean.  In  order  of  decreasing  albedo  they  are  (1)  dry,  cold 
snow,  wind  packed,  deep  drift,  p = 0.4  gf cm3,  (2)  5-crn  wet  new 
snow  over  multiyear  white  ice,  and  (3)  old  melting  snow,  28  cm 
thick.  Circles  are  calculated  albedo  of  semi-infinite  pure  snow  for 
diffuse  illumination,  with  grain  radii  to  match  observations  atX  = 0.9 
pm.  In  order  of  decreasing  albedo  they  are  (1)  r = 110  pm,  (2)  r = 
300  pm,  and  (3)  r = 1300  pm.  Plus  signs  are  for  snow  containing  the 
specified  concentrations  of  soot,  using  for  the  imaginary  index  of 
refraction  mim  (soot)  = 0.5  independent  of  wavelength;  = 1.8. 
(Taken  from  Figure  6b  of  WWII  but  with  soot  amounts  corrected 
according  to  footnote  3 of  WWII). 
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Fig.  7.  Effect  of  Mount  St.  Helens  ash  on  snow  albedo  for 
diffuse  incidence:  (a)  snow  grain  size  r - 100  pm:  (6)  r = 1000  pm. 
Ash  parameters  are  described  in  the  text. 


giving  no  color  to  the  snow;  but  in  order  to  mimic  a given 
concentration  of  soot,  the  andesite  concentration  must  be 
200  times  higher. 

For  this  review  article  we  have  done  a special  calculation 
of  the  effect  on  snow  albedo  of  Mount  St.  Helens  ash  (Figure 
7),  using  the  model  of  WWII,  for  diffuse  incident  radiation. 
The  ash  particle  size  distribution  is  taken  to  be  that  mea- 
sured from  aircraft  at  3200-m  elevation,  130  km  downwind  of 
the  volcano  on  May  19,  1980  [Hobbs  et  al.,  1981,  Figure  2c]. 
This  size  distribution  has  about  the  same  effect  on  snow 
albedo  as  does  a uniform  size  of  r = 5 pm.  The  real  refractive 
index  is  taken  as  that  of  andesite  (mre  = 1 .47)  from  Pollack  et 
al.  [1973].  The  imaginary  index  mim  (X)  was  measured  by 
Patterson  [1981]  for  0.3  £ X ^ 0.7  pm.  We  take  the  values  he 
gives  for  the  ash  which  feli  at  Bozeman,  Montana,  and  then 
interpolate  frem  his  value  at  0.7  pm  to  the  value  for  andesite 
at  1.2  pm.  For  X > 1.2  pm  we  use  the  andesite  values 
reported  by  Pollack  et  al.  In  the  visible  wavelengths  the 
Mount  St.  Helens  ash  is  measured  to  be  more  absorptive 
{mtm  ~ 4 x 10-3)  than  was  the  andesite  sample  of  Pollack  et 
al.  (mim  « 1 x 10~3)  and  also  to  be  slightly  reddish  colored 
rather  than  perfectly  grey. 

The  snow  albedo  in  Figure  7 is  reduced  for  X < 1 pm  by 
addition  of  ash.  At  longer  wavelengths  the  albedo  is  unaf- 
fected unless  the  ash  content  exceeds  0.1%,  in  which  case 
the  albedo  is  increased  for  X > 1.15  pm. 

The  eruption  of  Mount  St.  Helens  has  provided  an  oppor- 
tunity to  study  the  effects  of  volcanic  ash  on  snowmelt.  A 
curious  puzzle  reported  by  R.  Armstrong  (personal  commu- 
nication, 1980)  is  that  although  the  melting  rate  on  Mount 
Olympus  in  summer  1980  was  enhanced  by  the  Mount  St. 
Helens  ash,  the  formation  of  suncups  [Post  and  LaChapelle, 
1971,  pp.  71-73]  was  dramatically  inhibited.  The  solution  to 
this  puzzle  may  lead  to  further  insight  into  the  nature  of  the 
snowmelt  process. 
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Fig.  8.  Measurements  of  the  downward  flux  of  solar  radiation  in 
the  snow  at  Plateau  Station,  Antarctica.  The  horizontal  axis  is 
geometric  depth,  not  liquid-equivalent  depth.  (Figure  1 of  Schwerdu 
feger  and  Weller  [1977].) 


I.  Transmission  of  Light  Through  Snow 
1.  Spectrally  Integrated  Flux  Extinction 

Figure  8,  taken  from  Schwerdtfeger  and  Weller  [1977] 
shows  the  extinction  of  all-wave  solar  radiation  flux  in  clean, 
dry  snow  ( p , = 0.3  g cm-3)  of  the  Antarctic  Plateau.  The 
extinction  is  due  to  both  scattering  and  absorption.  It  does 
not  follow  an  exponential  decay  in  the  top  40  cm  or  so 
because  the  spectral  composition  of  the  downflux  changes 
rapidly  in  this  region.  Only  the  visible  radiation  penetrates 
deeply.  At  a depth  of  100  cm  (30  cm  liquid  equivalent),  1%  of 
the  incident  downflux  remains  unextinguished.  This  \%  is 
probably  entirely  blue  light,  concentrated  in  a narrow  wave- 
length band  about  X = 0.46  pm,  where  k,-  reaches  its 
minimum. 

The  absorption  of  solar  radiation  is  greatest  near  the 
surface  and  is  due  almcst  entirely  to  near-IR  radiation 
[Wiscombe  i.,id  Warren,  19806,  Figure  3]  because  most  of 
the  visible  light  eventually  reemerges  after  multiple  scatter- 
ings. Choudhury  [1981,  Figure  14]  has  calculated  the  solar 
flux  divergence  and  the  associated  heating  rates  in  the  snow 
using  an  atmospheric  radiation  model  coupled  to  his  snow 
albedo  model.  He  finds  the  heating  rate  at  the  surface  to  be 
about  20  times  that  at  5-nun  depth  for  r = 300  pm  but 
dependent  on  grain  size.  The  factor  is  larger  for  smaller 
grains  because  radiation  is  attenuated  more  rapidly  in  fine- 
grained snow. 

These  calculated  solar  heating  rates  will,  of  course,  be 
compensated  by  infrared  cooling  rates.  The  k,  for  thermal  IR 
radiation  is  much  larger  than  that  for  near-IR  radiation,  so 
the  cooling  is  more  concentrated  at  the  surface  than  is  the 
heating.  This  leads  to  the  observation  of  temperature  maxi- 
ma at  some  depth  below  the  surface  in  polar  snowfields 
during  the  sunlit  seasons.  The  temperature  maximum  at 
Pionerskaya  (Antarctica)  in  December  was  located  8 cm 
below  the  surface  [Schlatter,  1972,  Figure  2],  Similar  results 


are  found  in  Greenland  (W.  Bow,  personal  communication, 
1981).  The  depth  of  maximum  net  heating  rate  should  vary 
with  solar  zenith  angle,  approaching  the  surface  as  0o~*  90°. 

2.  Spectral  Flux  Extinction 

For  the  testing  of  theoretical  models,  monochromatic  flux 
extinction  measurements  are  more  useful  than  the  spectrally 
integrated  measurements  of  Figure  8 because  of  the  changing 
spectral  composition  with  depth  and  the  uncertain  spectral 
composition  of  the  incident  sunlight. 

Flux  extinction  measurements  are  more  difficult  to  do 
accurately  than  are  albedo  measurements  for  diffuse  inci- 
dence. (Aibedo  measurements  under  clear  sky  arc  also 
subject  to  error,  as  is  explained  in  section  J below.)  The 
presence  of  the  radiometer  in  the  snowpack  disturbs  the 
radiation  field  more  than  it  does  above  the  snowpack;  this 
disturbance  was  analyzed  by  Giddings  and  LaChapelle 
[1961].  A second  difficulty  which  is  more  severe  for  extinc- 
tion measurements  than  for  albedo  is  ‘leakage,’  the  detection 
of  unwanted  light  from  the  wings  of  the  filter  function  (T.  C. 
Grenfell,  personal  communication,  1979).  Leakage  is  least 
important  at  the  wavelength  of  smallest  absorption,  X = 0.47 
pm.  For  example,  one  may  nominally  be  measuring  radia- 
tion flux  at  X,  because  the  filter  function  is  centered  at  X|, 
but  some  smali  fraction  of  the  light  at  X2  in  the  wing  of  the 
filter  function  also  enters  the  detector.  Unfortunately,  the 
downflux  at  X2  deep  in  the  snowpack  may  be  orders  of 
magnitude  larger  than  that  at  X|,  so  that  most  of  the  light  one 
reports  at  X|  was  actually  light  of  the  wrong  wavelength. 
This  is  a problem  for  extinction  measurements  but  not  for 
albedo  measurements:  for  example,  whereas  the  ratio  of  blue 
albedo  (X  = 0.46  pm)  to  red  albedo  (X  = 0.7  pm)  may  be 
close  to  1.0  (see  Figure  1),  the  ratio  of  blue  downflux  to  red 
downflux  (for  r = 1 10  pm  in  Figure  9a)  may  be  —60  at  a 
depth  of  10  cm  (liquid  equivalent)  anu  —4000  at  20  cm. 

Figure  9 compares  measurements  with  calculations  of  ks 
(X).  Liljequist  [1956]  used  four  filters  to  study  the  visible 
spectrum  .•  i Antarctic  snow  whose  grain  size  he  measured  as 
r = 150  pm.  Kuhn  and  Siogas’  [1978]  (hereafter  KS) 
measurements  in  south  polar  snow  probably  had  good  spec- 
tral resolution  throughout  the  range,  whereas  Grenfell  and 
Maykufs  [1977]  (hereafter  GM)  spectral  resolution  became 
poorer  with  increasing  wavelength.  Deconvolution  was  used 
for  the  longer  wavelengths  in  an  attempt  to  improve  the 
resolution  (T.  C.  Grenfell,  personal  communication,  1981). 
The  data  of  GM  in  Figure  9 are  probably  reliable  for  X < 0.6 
pm. 

Albedo  had  also  been  measured  for  these  snow  samples, 
and  grain  size  was  obtained  by  WWI  and  WWII  by  matching 
the  calculated  and  observed  albedo  at  X a 0.9  pm,  where 
impurities  have  negligible  effect  on  albedo.  (This  was  done 
for  the  albedos  of  KS  and  GM;  Liljequist’s  albedos  indicate 
that  the  snow  at  Maudheim  was  uncontaminated.)  The  grain 
sizes  thus  obtained  (Figure  6)  were  r = 1300  pm  for  GM’s  old 
melting  snow  and  r =»•  1 10  pm  for  both  GM’s  ‘dry  compact 
snow’  and  KS’s  south  polar  snow  (r  = 110  pm  also  tolerably 
matches  Liljequist’s  albedos). 

The  model  parameters  necessary  to  explain  the  albedo 
agreed  with  the  parameters  necessary  to  eAplain  ks  only  for 
Liljequist's  measurements.  (However,  his  measurements  are 
unreliable  for  X a 0.6  pm  because  of  leakage  as  described 
above.)  The  presence  of  0.5  ppmw  of  soot  was  found 
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necessary  by  V/W|l  (after  taking  into  account  footnotes  3 
ami  4 ol  WWII)  to  explain  the  KS  albedo  measurements,  yet 
J iguie  9ti  sliows  that  KS  s flux  extinction  measurements  arc 
characteristic  of  pure  snow,  as  can  also  be  seen  in  Figure  1 1 
of  Choudhury  [1 981). 

I lie  calculations  (circles  and  plus  signs)  in  Figure  9 
correspond  to  the  circles  and  plus  signs  in  Figure  6;  the  plus 
signs  give  k,  for  the  soot  amount  necessary  to  explain  the 
albedos  of  GM.  Considering  only  X < 0.6  /cm  for  reasons 
mentioned  above,  the  measured  ks  falls  between  the  pure- 
snow  k,  and  the  sooty-snow  x,.  This  might  be  explained  if 
the  soot  concentration  at  the  surface  at  T-3  had  been  higher 
than  its  subsurface  concentration.  This  is  possible  if  the  soot 
was  the  result  of  combustion  of  heating  fuel,  because  the 
nearby  building  had  been  unoccupied  until  shortly  before  the 
measurements  were  made  (T.  C.  Grenfell,  persona!  commu- 
nication, 1981).  However,  this  argument  probably  cannot 
explain  the  data  of  Kuhn  and  Siogas;  it  is  difficult  to  see  how 
soot  could  be  concentrated  at  the  surface  1 km  from  the 
South  Pole  Station,  rather  than  uniformly  distributed. 


Fig.  9.  Asymptotic  flux  extinction  coefficient  for  snow  as  a 
function  of  wavelength  in  the  visible  spectrum.  Solid  lines  are 
measurements  of  Grenfell  and  Maykut  [1977,  Figure  3].  Squares  are 
Antarctic  measurements  of  Uljequlst  (1936)  for  snow  of  measured 
grain  size  r =»  130  /xm.  Bandwidfns  of  LUjequist's  filters  are  obtained 
from  his  Figure  48  and  plotted ' sre  as  the  positions  of  half-maximum 
transmittance.  Triangles  measurements  of  Kuhn  and  Siogas 
)I978|  at  South  Pole  Station.  (All  depth  values  Have  been  converted 
here  from  snow  depth  to  liquid-equivalent  depth.)  Circles  are 
calculations  using  WWI  model,  with  grain  radii  chosen  to  match 
albedo  at  X « 0.9  pm  for  snowpacks  shown  in  Figure  1 of  Grenfell 
and  Maykut  [1977].  These  are  the  same  model  snowpacks  described 
by  (he  top  and  bottom  lines  of  the  circles  In  Figure  6 of  this  paper. 
Plus  signs  are  model  calculations  for  snowpacks  corresponding  to 
the  top  and  bottom  lines  of  the  plus  signs  in  Figure  6. 
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Fig.  10  Dircct-bcam  albedo  of  a semi-infinite  snowpack  versus 
wavelength  for  several  values  of  direct-beam  zenith  angk  0 - cos' 1 
/i<>-  (Figerc  1 1 a of  WWI.  Reprinted  by  permis.ion  of  the  American 
Meteorological  Society.) 

These  discrepancies  indicate  the  need  for  more  spectrally 
detailed  measurements  of  both  k,(X)  and  a/X)  on  the  same 
snowpack.  As  Choudhury  [1981]  emphasized,  a good  theo- 
retical model  must  explain  both  albedo  and  extinction  mea- 
surements. With  the  currently  available  data  we  place  more 
emphasis  on  explaining  the  albedo  measurements  than  on 
explaining  the  k,  measurements,  because  the  latter  are  more 
susceptible  to  experimental  error, 

J.  Dependenll  of  Snow  Albedo  on  Sun  Ancle 
1 . Spectral  Snow  Albedo 

As  is  true  for  most  surfaces,  the  albedo  of  snow  increases 
as  the  sun  nears  the  horizon.  Figure  11a  of  WWI  is  repro- 
duced here  as  Figure  10,  showing  model-calculated  albedo 
for  four  zenith  angles.  The  albedo  is  predicted  to  be  sensitive 
to  zenith  angle  in  the  near  IR  but  not  in  the  visible,  in 
qualitative  agreement  with  measurements  of  Bryazgin  and 
Koptev  [1969],  However,  addition  of  trace  amounts  of 
impurities  to  the  snow  would  allow  a zenith  angle  depen- 
dence in  the  visible  as  well.  (These  calculations  arc  for  the 
extreme  case  of  100%  direct-beam  radiation.  Even  under 
dear  sky  the  diffuse  sky  radiation  would  make  the  observed 
zenith  angle  dependence  weaker  than  that  shown  in  Figure 
10.) 

The  reason  that  the  albedo  is  higher  for  low  sun  is  that  a 
photon  on  average  undergoes  its  first  scattering  event  closer 
to  the  surface  if  it  entered  the  snow  at  a grazing  angle.  If  the 
scattering  event  sends  it  in  an  upward  direction,  its  chance  of 
escaping  the  snowpack  without  being  absorbed  is  greater 
than  it  would  be  if  it  were  scab'ertd  from  deeper  in  the  pack. 
This  would  be  observed  even  if  tr  e ice  particles  scattered 
light  equally  in  all  directions.  Bui  the  phenomenon  is  greatly 
enhanced  by  the  extreme  asymmetry  of  the  scattering, 
whereby  scattering  within  a few  degrees  of  the  forward 
direction  is  much  more  probable  tnan  scattering  to  other 
angles. 

In  all  radiative  transfer  problems  for  a plane-parallel  slib, 
as  the  direction  of  incidence  goes  toward  grazing,  the  albedo 
becomes  increasingly  dominated  by  single  scattering.  As  the 
sun  goes  down,  the  shoulder  of  the  forward  peak  of  the 
single  scattering  phase  function  begins  to  emerge  more  and 
more  from  the  snowpack  (Figure  1 1),  whereas  when  the  sun 
was  higher  It  was  buried  because  of  its  5*-l0*  width.  Natural- 
ly. it  is  the  grains  at  the  surface  that  are  doing  most  or  the 
scattering  in  this  case.  Figure  1 1 is  purely  schematic,  be- 


Fig.  ! 1 . Polar  diagram  of  the  scattering  phase  function  of  a snow 
particle  at  the  surface  of  the  snowpack  for  three  solar  zenith  angles. 
This  shows  the  probability  that  a scattered  photon  will  go  into  any 
particular  direction.  For  display  purposes  this  phase  function  was 
calculated  for  an  unrealistically  small  snow  grain  size  (r  = 10  pm,  X 
= 5 p m).  The  asymmetry  of  the  phase  function  becomes  much  more 
extreme  as  r increases  or  X decreases. 

cause  it  was  calculated  for  spherical  ice  particles.  For 
hexagonal  columns  or  plates  the  22°  halo  may  contain  as 
much  energy  as  the  forward  peak  [Wendling  et  al.,  1979, 
Figure  S],  and  this  halo  begins  to  emerge  from  the  snowpack 
already  at  fib  “ 65°. 

With  very  low  sun  we  approach  a pure  single-scattering 
situation,  so  that  the  bidirectional  reflectance  becomes  just 
the  single-scattering  phase  function. 

At  the  same  time  the  assumption  of  sphericity  of  snow 
grains  breaks  down.  The  effects  of  angular  details  are 
smeared  out  with  multiple  scattering,  so  that  the  grains  can 
be  treated  as  spheres.  But  the  angular  effects  become 
dominant  a s -» 90°. 

Besides  the  breakdown  of  the  sphericity  assumptions  as  fib 
—»  90°,  the  delta-Eddington  approximation  also  becomes 
poor,  for  reasons  which  are  given  by  Wiscombe  [1977].  The 
disagreement  between  model  and  observation  as  -*  90°  is 
attributed  (below)  to  this  breakdown  of  delta-Eddington 
approximation. 

All  of  this  can  be  discussed  in  terms  of  standard  single- 
scattering and  multiple-scattering  theory.  There  is  no  need  to 
invoke  a separate  'surface  reflection.’  Specular  reflection  is 
an  acceptable  natural  part  of  the  single-scattering  pattern  of 
hexagonal  plates  and  prisms  and  other  ice  crystal  forms. 
Hence  if  one  would  use  a fully  correct  single-scattering 
phase  function  instead  of  the  spherical  Mie  phase  function, 
and  an  exact  radiative  transfer  method  for  multiple  scatter- 
ing, one  could  properly  model  the  albedo  as  6b  -*  90°. 

There  are  no  nomochromatic  measurements  against  which 
to  test  the  model.  This  is  unfortunate,  because  the  discrepan- 
cies in  spectrally  integrated  albedo  between  model  and 
observation  described  below  might  be  pinned  down  if  mea- 
surements would  be  made  at  discrete  wavelengths. 

2,  Spectrally  Integrated  Snow  Albedo  a {po) 

This  is  a controversial  subject  at  present.  The  reasons  for 
the  contusion  are  (1)  modeling  error,  (2)  instrument  error,  (3) 
inadequate  observation  of  snow  grain  size,  impurity  content, 
and  surface  roughness,  all  of  which  should  affect  the  slope  of 
a(p o),  and  (4)  inadequate  knowledge  of  the  spectral  distribu- 
tion of  the  incident  radiation. 

A number  of  clear-sky  measurements  of  spectrally  inte- 
grated snow  albedo  as  a function  of  the  cosine  Pq  of  the  solar 
zenith  angle  6b  are  plotted  in  Figure  12,  along  with  model 
results.  (The  only  reason  for  plotting  the  data  on  the  two 


separate  frames  (Figures  12a  and  12 b)  is  to  avoid  clutter.) 
Considering  first  Figure  12a,  we  note  that  the  model  calcula- 
tions of  Wiscombe  and  Warren  [1980b]  (for  the  atmospheric 
conditions  of  the  Antarctic  Plateau  and  r = 100  pm)  show  a 
much  weaker  dependence  of  a on  than  do  the  observa- 
tions by  Hubley  [1955]  on  the  Juneau  Icefield,  Korff  et  al. 
[1974]  in  Colorado,  and  Rusin  [1961]  in  Antarctica.  (Note 
that  the  model  does  not  account  for  the  sphericity  of  the 
atmosphere,  which  probably  makes  it  unreliable  for  /zq  <0.1 
or  so.) 

The  responses  of  all  commercial  radiometers  deviate  from 
a proper  ‘cosine  law’  ] Liljequist , 1956;  Dirmhirn  and  Eaton, 
1975].  They  are  usually  less  sensitive  at  large  incident  zenith 
angles.  If  not  corrected  for,  this  causes  albedos  at  low  sun  to 
be  overestimated.  This  is  because  the  reflected  radiation  is 
more  diffuse  than  the  incident  radiation.  The  steep  depen- 
dence of  a on  po  reported  by  Rusin  and  by  Hubley  might  be 
explained  as  this  type  of  error.  But  the  measurements  of 
Korff  et  al.  must  be  considered  reliable  because  they  cali- 
brated their  instrument  and  applied  a correction  at  large 
zenith  angles. 

Hubley’s  results  show  an  interesting  hysteresis,  with 
higher  albedo  in  the  morning  than  afternoon  at  the  same 
zenith  angle.  Hubley  speculated  that  this  might  be  attributed 
to  specular  reflection  from  fimspiege!  in  the  morning  before 
the  melting  began. 

In  Figure  12b  are  plotted  several  Antarctic  measurements 
and  calculations.  Carroll  and  Fitch' s [1981]  (hereafter  CF) 
measurements  extended  to  larger  zenith  angles  than  anyone 
else  has  reported  and  show  a much  steeper  dependence  of  d 
on  po  than  Li(jequis:  [1956]  found.  Liljequist’s  instruments 
had  been  extensively  calibrated  for  the  dependence  of  their 
response  on  both  the  zenith  and  the  azimuth  solar  angles 
[Liljequist,  1956,  pp.  45-55].  CF  did  not  report  such  calibra- 
tions for  their  instrument,  so  it  is  possible  that  their  dramatic 
disagreement  with  Liljequist  is  at  least  partly  caused  by  an 
experimenta’  bias  in  CF’s  measurements. 

CF’s  measurements  were  complicated  by  surface  irregu- 
larity. The  faces  of  oriented  sastrogi  present  different  angles 
to  the  sun  as  it  moves  around  the  horizon.  Accordingly,  CF 
noted  a diurnal  cycle  of  albedo,  as  had  Kuhn  a, id  Siogas 
[1978],  with  higher  values  when  the  sun  was  oriented  parallel 
to  the  sastrugi,  so  that  the  effective  zenith  angle  was  larger. 
The  measurements  of  CF  plotted  in  Figure  12b  are  the 
averages  of  four  solar  azimuths  6 hours  apart.  A repre- 
sentative error  bar  is  drawn  on  one  of  the  points. 

There  are  actually  several  effects  of  surface  roughness  on 
albedo.  In  a field  of  randomly  oriented  surface  roughness 
features  (or  in  the  daily  average  at  the  south  pole)  when  the 
sun  is  low,  the  effective  zenith  angle  is  always  smaller  for  a 
rough  surface  than  for  a flat  surface  (section  5c  of  WWI). 
This  means  that  CF’s  observed  dependence  of  a on  po  would 
be  even  steeper  if  p<>  were  taken  as  the  effective  value 
instead  of  the  flat  surface  value,  and  their  results  would 
deviate  even  more  dramatically  from  the  results  of  other 
experimenters  plotted  in  Figure  12. 

However,  this  is  further  complicated  by  the  fact  that  the 
effective  zenith  angle  is  a function  of  wavelength.  In  order 
for  surface  roughness  to  moderate  the  zenith  angle  effect  its 
typical  length  and  depth  must  not  be  much  smaller  than  the 
average  penetration  depth  of  light  into  snow.  Surface  irregu- 
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Fig.  12.  Spectrally  integrated  clear-sky  snow  albedo  as  a function  of  r0  according  to  various  investigators.  The 
resu.t  of  the  Wiscombe  and  Warren  [19806]  Antarctic  model  is  plotted  in  both  frames  for  reference.  The  lower  two  plots 
in  Figure  126  are  for  Antarctic  planetary  albedo. 


larities  on  smaller  scales  are  less  important  because  a 
significant  amount  of  radiation  penetrates  right  through  a 
peak  and  emerges  on  the  other  side.  Thus  whereas  surface 
roughness  of  millimeter  seals  affects  the  reflection  of  visible 
light  by  soil  particles,  it  does  not  affect  the  reflection  of 
visible  light  by  snow. 

Besides  altering  the  effective  zenith  angle,  surface  rough- 
ness will  also  decrease  the  albedo  further  by  trapping  some 
radiation  in  the  troughs,  as  described  in  section  E5 b. 

In  addition  to  these  systematic  effects  there  is  a sampling 
error  which  arises  in  albedo  measurements  of  a rough 
surface.  One  should  ideally  make  measurements  from  high 
towers  [Langleben,  1968]  so  as  to  get  a representative  view 
of  the  surface.  Close  to  the  ground  the  field  of  view  may 
contain  only  one  or  a few  sastrugi  and  will  thus  contain  more 
or  less  than  its  fair  share  of  shadow,  depending  on  solar 
azimuth. 

The  delta-Eddington  snow  albedo  model  of  WWI  was 
coupled  to  the  atmospheric  radiation  model  of  Wiscombe  for 
cloud-free  January  conditions  of  the  Antarctic  Plateau  [Wis- 
combe  and  Warren,  19806]  The  zenith  angle  dependences  of 
both  surface  albedo  and  planetary  albedo  are  shown  in 
Figure  12 b.  The  model  adequately  explains  the  a(po)  ob- 
served by  LUjeqnist  [1956]  but  clearly  fails  to  reproduce  that 
of  CF.  (The  reduction  of  December-January  albedo  at  large 
zenith  angles  which  shows  up  in  Ltljequist’s  daily  averages 
was  only  observed  in  the  afternoon  and  was  attributed  by 
Liljequist  to  metamorphism  rather  than  to  zenith  angle.  This 
is  a possible  complicating  factor  in  all  of  these  studies.) 

Choudhury  and  Chang  [1981]  added  'surface  reflection’  to 
their  model,  as  described  in  section  E4  above,  in  order  better 
to  match  the  November  measurements  of  Liljequist.  But  in 
order  to  get  even  this  slight  steepening  of  the  a(po)  shown  in 
Figure  126  they  had  to  use  what  seems  an  unreasonably 
small  slope  variance  s1  = 0.01,  which  would  be  more 
appropriate  for  a glazed  crust  than  for  the  dry  snow  grains  of 
Antarctica.  The  point  we  wish  to  emphasize  here  is  that  the 
addition  of  even  this  extreme  amount  of  surface  reflection  is 


dramatically  inadequate  to  bring  the  delta-Eddington  model 
in  line  with  CF’s  observations.  If  these  observations  must  be 
matched  instead  of  those  of  Liljequist,  a problem  with  the 
delta-Eddington  model  is  indicated  that  cannot  be  corrected 
by  adding  an  ad  hoc  surface  reflection. 

Delta-Eddington  is  known  to  underestimate  the  albedo  of 
a plane-parallel  layer  at  large  zenith  angles.  It  agrees  very 
well  with  the  exact  doubling  calculations  [see  Joseph  et  al.. 
1976,  Figure  3]  for  po  a 0.4  (So  £ 66°),  but  albedo  errors  of 
up  to  10%  can  occur  at  po  - 0.1  (fib  = 84°)  for  particular 
values  of  g and  w.  In  order  to  quantify  the  errors  of  the  delta- 
Eddington  method  for  snow  at  low  sun,  it  will  be  necessary 
to  do  some  albedo  calculations  which  give  a more  exact 
account  of  the  radiative,  transfer  process. 

3.  Spectrally  Integrated  Planetary  Albedo 

Model  calculations  of  spectrally  integrated  planetary  albe- 
do over  the  Antarctic  Plateau  plotted  in  the  lower  part  of 
Figure  126  show  planetary  albedo  to  increase  slightly  from 
Po  = 0.6  to  fio  = 0.3  because  of  the  p$  dependence  of  surface 
albedo.  But  when  the  sun  goes  even  lower,  a second  effect 
comes  into  play  to  reduce  the  albedo.  This  is  the  increased 
H20  and  O3  absorption  in  the  atmosphere  due  to  the  longer 
slant  path  of  the  sunlight. 

The  zenith  angle  dependence  of  planetary  albedo  of  the 
snow-atmosphere  system  is  being  measured  by  Stowe  et  al. 
[1980,  personal  communication,  1980]  using  the  scanning 
radiometer  of  the  Nimbus  7 satellite.  For  that  work  all  snow- 
covered  scenes  at  the  same  zenith  angle  are  considered 
equivalent,  regardtes  of  geographic  location.  However,  for 
the  1 month  of  data  available  so  far  (November  1978),  nearly 
all  of  the  scenes  were  over  Antarctica.  These  observations  of 
Stowe  et  al.  plotted  as  the  histogram  in  Figure  126  show  a 
steeper  decrease  in  albedo  as  po  decreases.  It  is  probably 
premature  to  comment  on  the  differences  between  Stowe’s 
observations  and  the  model  predictions,  since  many  factors 
remain  to  be  considered,  but  it  is  noteworthy  that  both 
model  and  observation  agree  on  a range  of  about  70-73%  for 
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Fig.  13.  Effect  of  cloud  cover  on  spectrally  integrated  snow 
albedo.  The  difference  between  observed  and  calculated  spectrally 
integrated  albedos  is  plotted  versus  incident  irradiance  (which  is 
here  approximately  proportional  to  atmospheric  transmissivity, 
because  the  measurements  were  all  made  between  1100  and  1400) 
for  snow  in  the  Cascade  Mountains.  'Calculated'  albedos  were 
calculated  by  Grenfell  el  at.  {19811  using  observed  spectral  albedos 
with  a fixed  incident  spectral  distribution  characteristic  of  an 
intermediate  cloud  thickness.  The  straight  line  is  a least  squares  fit 
to  the  data  points.  (Figure  5 of  Grenfell  el  at.  [1981].  Reprinted  by 
permission  of  Elsevier.) 


pio  - 0.2.  Of  all  the  surface  types  examined  by  Stowe  et  al., 
snow  was  the  only  one  in  which  the  atmospheric  path  length 
effect  was  able  to  dominate  over  the  surface  albedo  effect.  In 
order  to  be  consistent  with  these  planetary  albedos,  the  po 
dependence  of  the  snow  surface  albedo  would  have  to  be 
even  smaller  than  the  modest  dependence  given  by  the  delta- 
Eddington  model. 

In  summary,  (1)  the  delta-Eddington  method  underes- 
timates albedo  at  low  sun  (the  magnitude  of  this  underesti- 
mate could  be  determined  by  doing  a more  exact  calcula- 
tion), and  (2)  the  Antarctic  planetary  albedo  measurements 
of  Stowe  et  al.  and  the  Antarctic  surface  albedo  measure- 
ments of  Carroll  and  Fitch  are  mutually  inconsis*..  it. 

Our  detailed  discussion  here  of  extreme  zenith  angles  (Bo 
~ * 90°)  is  important  for  focusing  on  discrepancies  between 
model  and  observation,  but  it  has  very  limited  relevance  for 
the  snow  energy  budget.  Even  at  the  north  and  south  poles, 
only  3%  of  the  annual  irradiance  is  received  at  po  < 0. 1 ($>  > 
84.3°),  and  only  13%  at  p<>  < 0-2  (0o  > 78.5°).  At  other 
latitudes  these  percentages  are  even  smaller. 

K.  Effects  of  Cloud  Cover  on  Snow  Albedo 

1.  Monochromatic  Albedo 

The  only  effect  of  cloud  cover  on  monochromatic  albedo  is 
to  diffuse  the  radiation,  changing  the  effective  zenith  angle. 
The  effective  zenith  angle  for  purely  diffuse  radiation  is 
about  50°.  Thus  as  shown  in  Figure  12  of  WWI,  the  interposi- 
tion of  a cloud  layer  between  sun  and  snow  causes  spectral 
snow  albedo  to  increase  for  fib  < 50°  and  to  decrease  for  6b  > 
50%  the  latter  being  the  normal  situation  over  snow-covered 
surfaces. 

2.  Spectrally  Integrated  Albedo 

Cloud  cover  is  normally  observed  to  cause  an  increase  in 
spectrally  integrated  snow  albedo.  This  is  because  clouds 
absorb  the  same  near-infrared  radiation  that  snow  would 
absorb,  leaving  the  shorter  wavelengths  (for  which  snow 
albedo  is  higher  (cf.  Figure  1))  to  penetrate  to  the  surface. 


This  was  pointed  out  by  Liljequist  [1956,  p.  88],  Grenfell  and 
Maykut  [1977,  p.  457],  and  WWI  (section  5 d).  This  'spectral 
shift’  effect  (named  by  CF)  dominates  the  contrary  ’zenith 
angle  alteration’  effect  discussed  above  but  is  moderated  by 
it,  as  seen  in  Table  173  of  Rusin  [1961]. 

An  example  of  the  spectral  shift  effect  is  shown  in  Figure 
13,  which  is  taken  from  Grenfell  et  al.  [1981].  The  difference 
dots  ~ dcaic  is  plotted,  where  da ic  was  based  on  the  spectral 
albedo  measurements  and  was  calculated  by  Grenfell  et  al. 
assuming  a fixed  spectral  distribution  of  incident  radiation 
characteristic  of  an  intermediate  cloud  thickness.  Under  a 
thick  overcast  the  albedo  is  found  to  be  1-9%  higher  than  it 
would  have  been  under  an  intermediate  cloud  condition. 

An  interesting  exception  to  the  general  rule  that  snow 
albedo  is  higher  under  cloud  cover  has  been  reported  by  CF 
to  occur  at  very  low  sun  and  can  be  directly  traced  to  the 
much  steeper  dependence  of  a on  po  that  they  found  at  these 
sun  angles  (Figure  12 b).  At  'he  low  solar  elevations  encoun- 
terd  by  CF  at  the  south  po..,  the  change  in  effective  zenith 
angle  caused  by  clouds  is  apparently  more  than  able  to 
compensate  for  the  spectral  shift  effect. 

L.  Bidirectional  Reflectance  of  Snow 


The  radiation  reflected  by  a snow  surface  is  not  distributed 
uniformly  into  al!  angles.  If  it  were,  then  the  principle  of 
reciprocity  [Siegel  and  Howell.  1972]  would  indicate  that 
snow  albedo  was  independent  of  6b-  The  pattern  of  reflection 
is  described  by  the  anisotropic  reflectance  function/,  which 
depends  on  the  source  zenith  angle  $o,  detector  nadir  angle 
ff,  and  the  relative  azimuth  <j>'  - <fo.  The  bidirectional 
reflectance  R is  the  product  of  the  albedo  and/,  as  given  by 
equation  (1).  If  the  surface  conditions  depend  or.  azimuth,  as 
is  the  case  for  the  oriented  sastrugi  at  the  south  pole,  / is 
actually  a function  of  the  two  individual  azimuth  angles  <k> 
and  <f>'  rather  than  merely  of  their  difference,  but  this 
complication  is  usually  ignored;/may  also  depend  on  snow 
grain  size  and  snow  depth. 

Satellite  detectors  with  narrow  fields  of  view  measure  R at 
only  one  or  a few  angles.  In  order  to  obtain  the  planetary 
albedo  from  individual  satellite  measurements  it  is  thus 
necessary  to  have  prior  knowledge  of/.  Measurements  of/ 
have  been  made  at  the  surface  by  Middleton  and  Mungall 
[1952]  and  Dirmhirn  and  Eaton  [1975],  from  low-flying 
aircraft  by  Griggs  and  Marggrafl  1967]  and  Salomonson  and 
Marlatt  [19686]  (with  more  details  by  Salomonson  and 
Marlatt  [1968a]),  and  from  satellite  by  Stowe  et  al.  [1980, 
personal  communication,  1980]. 

Because  the  single-scattering  phase  function  varies  with 
wavelength,/  could  additionally  be  a function  of  wavelength. 
However,  Griggs  and  Marggraf  [1967]  have  found /to  be 
independent  of  X,  at  least  for  0.44  ^ X £ 0.96  pm.  If  this 
holds  also  at  longer  wavelengths  (which  seems  unlikely),  it 
means  that  / can  be  obtained  from  spectrally  integrated 
measurements.  1 

Figure  14  shows  the  results  of  the  two  most  comprehen-  j 

sive  sets  of  measurements.  (The  other  reports  cited  above  | 

included  measurements  at  only  a few  angles.)  Dirmhirn  and  I 

Eaton  [1975]  measured /for  a melting  snowpack  five  times  | 

(i.e.,  at  five  sun  angles)  on  a single  afternoon  in  Utah.  (They  § 

reported  that  snow  appeared  glazed  at  sunrise  and  sunset,  { 


1 . Complete  Description  (for  Satellite 
Measurements) 
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Fig.  14.  Anisotropic  reflectance  functions  for  snow  at  the  surface  [from  Dirmhirn  and  Eaton,  1975,  Figure  3]  and  at 
the  top  of  the  atmosphere  (replotted  from  data  of  Stowe  et  at.  (1980,  personal  communication,  1980]). 


but  since  it  was  melting  in  the  afternoon  when  measurements 
were  taken,  we  doubt  that  glazing  affected  the  measure- 
ments.) A separate  diagram  is  given  for  each  solar  zenith 
angle  6b;  at  each  ./  was  measured  at  eight  azimuths  <j>'  - 
for  four  nadir  angles  O’:  0°,  30°,  60°,  and  90°.  In  each 
diagram, /is  proportional  to  the  distance  front  the  center  of 
the  inner  circle  to  the  appropriate  point  on  the  curves,  but 
for  use  it  must  be  normalized  for  each  diagram  so  that 
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flu 0.  m'.  P - <M  V-'  d/i'  dtp  = ! (4) 

hemisphere 

The  solar  beam  is  taken  to  be  incident  from  below.  There  is 
enhanced  reflectance  in  the  forward  direction,  and  the 
anisotropy  inc  eases  as >o  decreases  (i.e.,  as  the  sun  goes 
down). 

Stowe  et  al.  (1900]  have  taken  scanner  data  from  the 
Nimbus  7 satellite  and  grouped  them  into  bins  of  (6b,  O',  <f>'  - 
<fo).  Since  the  satellite  did  not  examine  a single  target  at  all 
angles,  the  measurements  from  all  snow  surfaces  on  earth 
were  treated  together  and  assumed  to  have  the  same  /. 
However,  at  the  time  of  this  pilot  experiment  (November 
1978)  most  of  the  earth’s  sunlit  snow  area  was  actually  in  the 
Antarctic,  so  the  results  of  Stowe  et  al.  (top  of  Figure  14)  are 
really  representative  of  a snow  surface  uncontaminated  by 
vegetation.  They  show  less  anisotropy  than  do  the  measure- 
ments at  the  surface.  This  difference  may  partly  be  due  to  the 
diffusing  effect  of  the  intervening  atmosphere,  but  there  may 
also  be  an  actual  difference  in /between  the  two  snow  types 
(dry,  clean,  fine-grained  Antarctic  snow  versus  old,  melting, 
somewhat  glazed  mid-latitude  snow).  These  results  were 
from  only  a small  sample  of  the  satellite  data  and  will  be 
augmented  in  the  future. 

The  bidirectional  reflectance  of  snow  has  not  been  mod- 
eled since  the  work  of  Barkstrom  and  Querfeld  [1975] 
reviewed  m section  E2  above.  The  doubling  [Hansen,  1969] 
or  discrete  ordinates  [Stamnes  and  Swanson,  1981]  radiative 


transfer  solutions  would  be  useful  in  this  connection,  al- 
though surface  roughness  may  be  as  important  as  doing  the 
plane-parallel  radiative  transfer  correctly,  because  of  its  two 
effects  discussed  above:  altering  the  effective  zenith  angle  by 
putting  some  of  the  snow  in  shadow  and  trapping  light  by 
multiple  reflections  among  troughs  and  peaks. 

2.  Azimuthally  Averaged  Bidirectional 
Reflectance  (for  Flux  Calculations) 

For  radiation  budget  calculations  which  produce  only 
fluxes  and  not  intensities,  the  full  BRDF  is  not  needed.  One 
only  requires  the  albedo  and  the  azimuthally  averaged 
anisotropic  reflectance  factor  r: 


flOo,  O' 


1 f2’ 

' > - ^ Jo  ' 


Mh  O’,  p - <M  dp 


S.  G.  Warren  and  W.  J.  Wiscombe  (unpublished  data,  1980) 
derived  a parameterization  for  / from  observations  of  snow 
reflectance.  They  then  used  it  to  apportion  the  model- 
calculated  snow  albedo  among  the  reflection  nadir  angles  O'. 

In  Figure  15  are  plotted  values  of  / they  obtained  from 
azimuthal  averaging  of  measurements  of/from  four  reports. 
Each  set  of  connected  points  corresponds  to  a particular 
zenith  angle  6b.  A straight  line  was  fit  to  each  of  these  10 
plots: 


/=  1 + b(ti'  - I) 


(5) 


(The  data  in  Figure  15  have  been  normalized  so  that  the  least 
squares  lines  all  pass  through  the  point  (1,  1).) 

The  slopes  b of  these  lines  are  plotted  versus  the  solar 
zenith  cosine  mt  in  Figure  16.  They  show  approximately  a 
linear  dependence  on  iiq,  so  the  resulting  parameterization 
for/is 


filMt,  ti')  = [3/(3  - />)][!  + b(ii'  - 1)] 


(6) 
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Fig.  15.  Azimuthally  averaged  anisotropic  reflectance  function/ 
of  snow  as  a function  of  nadir  angle  cosine  /z',  scaled  so  that  j (1)  = 
1.0  for  the  least  squares  linear  fit  of  each  data  set.  The  data  points 
were  obtained  by  azimuthally  averaging  charts  published  in  four 
reports,  as  follows:  Salomonson  and  Marlatt  (1968a,  b):  Measure- 
ments of  /for  snow  (depth  unspecified)  in  Utah  and  Wyoming  from 
aircraft  120  m above  the  surface,  wavelength  channel  0.55-0.85  /an. 
Results  were  very  similar  for  the  2-  to  4-pm  channel.  Middleton  and 
Mungall  (1952]:  Measurements  at  the  surface  using  artificial  incan- 
descent lighting.  Plotted  here  are  results  derived  from  their  Figure 
3a  for  wind-packed  snow.  The  anisotropy  exhibited  by  this  snow 
sample  was  greater  than  for  new  snow  fallen  in  calm.  Observations 
were  made  only  in  the  plane  of  incidence,  so  there  are  insufficient 
observations  to  do  an  azimuthal  average  for  any  zenith  angles 
except  0o  = 0°,  which  is  the  case  plotted  here.  Dirmhirn  and  Eaton 
(1975]:  Measurements  made  at  the  surface  under  natural  sunlight, 
April  afternoon  in  Utah,  melting  snow.  Griggs  and  Marggraf[l%7, 
p.  153]:  Measurements  from  aircraft,  120  m above  the  surface  of  a 
snow-covered  lake,  Oregon,  snow  depth  of  ■=!  m.  Measurements 
were  made  at  coarse  azimuthal  resolution  (A = 60°)  but  at  eight 
discrete  wavelengths  between  0.44  and  0.96  jun; /appeared  to  be 
independent  of  wavelength. 


where  h = 1.07/io  - 0.84.  The  prefix  3/(3  - b)  is  a 
normalizing  factor,  chosen  so  that  (4)  is  satisfied.  The 
dependences  of/on  p o and  p are  different,  so  (6)  does  not 
satisfy  the  reciprocity  principle  and  is  only  a first  attempt  to 
parameterize  available  data. 

The  measurements  used  to  obtain  (6)  were  made  under  a 
variety  of  snow  conditions  with  unspecified  grain  size,  and  it 
is  likely  that  (6)  is  not  really  correct  for  all  wavelengths  and 
grain  sizes.  However,  Wiscombe  and  Warren  [1980b]  found 
that  their  radiant  flux  calculations  were  actually  insensitive 
to  the  functional  form  of  /;  they  obtained  almost  identical 
fluxes  whether  they  used  (6)  or,  alternatively,  assumed 
isotropic  reflectance.  The  largest  differences  in  fluxes  at  the 
surface  occurred  in  the  visible  wavelengths  and  were  at  most 
0.6%. 


M.  Thermal  Infrared  Emission  From  Snow 
1.  Emissivity 

We  have  taken  many  pages  to  discuss  the  solar  reflectance 
of  snow,  yet  its  infrared  emission  deserves  only  a few 
paragraphs.  The  reason  for  this  is  that  the  infrared  emissivity 
of  snow  is  quite  insensitive  to  snowpack  parameters.  For 


many  purposes  one  can  simply  assume  a snow  emissivity  of 
about  99%,  as  was  measured  by  Griggs  [1968]. 

Berger  [1979]  adapted  the  derivations  of  Bohren  and 
Barkstrom  [1974]  to  the  limit  of  large  absorption  as  described 
in  section  E2  above.  His  calculated  emissivity  (Figure  17)  is 
independent  of  grain  size  but  decreases  slightly  with  increas- 
ing density.  This  calculated  density  dependence  may  depend 
on  the  particular  regular  array  of  spheres  assumed  by 
Berger,  as  noted  above. 

The  albedo  model  of  WWI  was  used  to  calculate  emissivi- 
ty (Figures  8b  and  lib  of  WWI)  by  virtue  of  KirchhofFs  law 
[Siegel  and  Howell,  1972,  p.  70]: 


40b,  X)  = 1 - aAOo,  X) 

where  as  is  the  albedo  for  zenith  angle  fib  and  e *s  the 
directional  thermal  emissivity  at  the  same  wavelength  into 
the  viewing  (nadir)  angle  0b • Figure  18  shows  the  calculated 
dependence  of  snow  emissivity  on  wavelength,  averaged 
over  emission  angle  for  five  grain  sizes  (Figure  18a)  and  for  a 
single  grain  size  r = 100  /zm  at  four  emission  angles  (Figure 
18b).  Since  these  calculations  were  done  without  a near-field 
correction,  we  do  not  show  results  beyond  40- run  wave- 
length, and  the  results  plotted  may  already  be  somewhat 
inaccurate  at  X = 20  /zm,  especially  for  r = 50  /zm. 

The  emissivity  is  sensitive  to  grain  size  only  at  certain 
wavelengths.  In  particular,  it  is  completely  insensitive  to 
grain  size  over  most  of  the  Planck  function  for  normal 
terrestrial  temperatures,  centered  near  X = 10  /zm.  The  Mie 
results  thus  support  Berger's  [1979]  assumption  of  no  trans- 
mission of  IR  radiation  through  ice  spheres  in  this  spectral 
range,  which  led  to  his  prediction  that  emissivity  would  be 
independent  of  grain  size. 

The  model  of  WWI  does  not  compute  a dependence  of  e 
on  />,.  However,  because  of  the  small  penetration  depth  of 
thermal  IR  radiation  in  ice,  ‘surface  roughness’  on  the  scale 
of  millimeters  may  affect  emissivity.  Berger's  [1979]  depen- 
dence of  e on  p,  is  actually  caused  only  by  the  relation  he 
obtained  between  ps  and  ‘projected  areas  of  surface  and 
subsurface  particles’  (his  equations  27  and  28),  which  seems 
actually  to  be  a measure  of  surface  roughness. 

Because  of  the  large  mm  of  ice  throughout  the  thermal 
infrared,  only  a very  thin  layer  of  snow  is  already  effectively 


Fig.  16.  Slopes  (b  in  equation  (5))  of  the  lines  in  Figure  15  are 
plotted  here  versus  Mo-  The  line  is  a least  squares  fit  to  the  points. 
Symbols  refer  to  the  dz*a  sources  given  in  Figure  15. 
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semi-infinite,  so  e is  independent  of  snow  layer  thickness. 
Small  amounts  of  impurities  will  also  not  affect  e. 

Figure  I8i>  shows  that  the  emissivity  is  near  100%  for 
overhead  viewing  (#<*,  = 0°)  but  is  significantly  lower  for 
large  nadir  angles.  Thus  the  satellite  viewing  angle  should  be 
taken  into  account  when  inferring  snow  temperatures  from 
satellite  infrared  channel  (\  - 11  pm)  emission  measure- 
ments. The  results  in  Figure  18  were  obtained  using  the 
delta-Eddington  approximation,  so  for  6m  = 80°  the  true 
emissivity  may  actually  be  lower  than  that  shown. 

2.  Brightness  Temperature 

If  e < 1,  the  brightness  temperature  TB  at  a single 
wavelength  will  be  smaller  than  the  true  temperature  T.  The 
difference  TB  - T depends  not  only  on  e but  also  on 
wavelength  because  of  the  variation  of  the  Planck  function 
with  wavelength.  J.  Dozier  (personal  communication,  1981) 
has  calculated  TB  - T for  3 s X £ 14  pm  for  snow  near  the 
melting  temperature,  plotted  in  Figure  19.  The  assumption 
that  snow  is  a blackbody  could  lead  to  an  underestimate  of 
temperature  by  as  much  as  2.5  K at  X = 14  pm  and  ff  = 75°. 
However,  the  estimation  of  snow  temperature  from  satellites 
will  likely  be  subject  to  more  error  from  uncertainty  in 
atmospheric  transmissivity  than  from  uncertainty  in  snow 
emissivity. 

N.  Remote  Sf.nsinc  of  Snow 
I.  Snowpack  Properties  From 
Albedo  Measurements 

Dozier  et  al.  (1981]  have  used  the  model  of  WWI  to 
calculate  snow  albedos  integrated  over  channels  1 and  2 of 
the  NOAA  Tiros  N satellite  (0.5-0.7  pm  and  0.7-1. 0 pm, 
respectively).  The  hope  is  to  deduce  grain  size  from  a near- 
IR  channel,  where  depth  and  contaminants  have  no  effect  on 
albedo,  and  then  use  the  deduced  grain  size  together  with  the 
channel  1 data  to  infer  snow  water  equivalent  depth  below 
some  threshold  value  around  100  mm.  Among  the  difficulties 
in  this  approach  are  (1)  the  conversion  of  bidirectional 
reflectance  to  albedo,  (2)  the  poor  location  of  channel  2 for 
this  purpose  (an  ideal  channel  would  be  located  in  the  region 
1.0-1. 2 pxa,  where  the  sensitivity  of  albedo  to  grain  size  is 
greatest),  and  (3)  the  tact  that  visible  albedo  reduction  can  be 
due  to  impurities  as  well  as  to  thinning  of  the  snowpack. 

Dozier  et  al.  were  apparently  able  to  detect  the  thinning  of 
the  snowpack  at  the  end  of  the  melting  season  on  some 
Canadian  lakes.  Figure  20  is  taken  from  their  work.  On  the 


3 5 7 D 20  30  40 

WAVELENGTH  (pm) 


ion  on 

I ittni 


Fig.  17.  Emissivity  as  a function  of  wavelength  for  various  snow 
densities,  according  to  the  model  of  Berger  [1979].  (Figure  7 of 
Berger  [1979].) 
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Fig.  18.  Thermal  infrared  emissivity  of  snow  as  a function  of 
wavelength,  according  to  the  model  of  WWI.  (a)  Hcmispherically 
averaged  emissivity  for  snow  grain  sizes  r = 50, 100,  200, 500,  1000 
pm.  ( b ) Directional  emissivity  for  four  detector  nadir  angles. 


earlier  date,  April  12,  the  snow-covered  lakes  are  brighter 
than  the  intervening  forest  in  both  channels.  Two  weeks 
later,  much  of  the  snow  had  melted,  so  that  on  April  27  there 
was  less  contrast  in  both  channels  between  lake  and  forest. 
On  both  of  these  dates  it  is  apparent  that  for  snow  the  near- 
IR  reflectance  is  smaller  than  the  visible  reflectance,  in 
agreement  with  Figure  1.  The  forest  exhibits  the  opposite 
behavior,  in  agreement  with  the  known  spectral  reflectance 
of  green  plants  [Carer,  1980].  Hopefully,  in  the  future  this 
work  will  be  combined  with  ground  truth  measurements  of 
snow  depth. 

Some  experiments  at  the  surface  which  are  relevant  to  this 
problem  are  being  done  by  O'Brien  and  Koh  [1981].  They 
observed  the  change  in  spectral  reflectance  (using  a few 
narrow-band  filters)  as  a thick  snow  cover  decayed,  docu- 
menting the  transition  from  the  spectral  reflectance  of  snow 
to  the  spectral  reflectance  of  grass.  As  is  expected  (Figure  13 
of  WWI),  the  underlying  surface  first  begins  to  ‘show 
through’  in  the  visible  wavelengths  and  is  evident  in  the  near 
IR  only  when  the  snow  cover  is  much  thinner.  Because  of 
the  somewhat  crude  experimental  setup,  these  results  are  as 
yet  only  qualitative. 

Detection  of  the  dust  content  of  snow  has  been  attempted 
by  Sydor  et  al.  [1979]  for  polluted  snow  in  Duluth  harbor. 
However,  dust  was  reported  as  rates  of  deposition  rather 
than  as  weight  fractions  in  snow,  so  no  direct  comparison 
with  a model  can  be  made. 

Matson  and  Wiesnet  [1981]  describe  the  routine  global 
mapping  of  snow  cover  using  visible-channel  imagery  from 
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Fig.  19.  Difference  between  brightness  temperatures  TB  and 
snow  temperature  T as  a function  of  wavelength  (a)  for  three 
different  snow  grain  sizes  at  viewing  angle  ff  = 60°  and  (6)  for  snow 
grain  radius  r = 500  /an  at  five  different  viewing  angles.  Figure  from 
J.  Dozier  (personal  communication,  1981). 

NOAA  satellites.  This  monitoring  program  has  been  useful 
in  describing  the  interannual  variability  of  regional  snow 
cover. 

2.  Snow  Cloud  Discriminator  {1.6  tun) 

Measurements  from  aircraft  of  spectral  reflected  intensity 
in  the  1.4-  1.8-/un  wavelength  region  have  shown  that  clouds 
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Fig.  20.  Traces  of  reflectance  measured  bom  satellite  by  NOAA 
saving  radiometer  across  three  lakes  in  Saskatchewan  on  two 
different  days  in  visible  channel  i (0.5-0.7  tun)  and  near-IR  channel 
2 (0.7-1.0  pm).  (Figure  10  of  Dozier  eta!.  (1981].) 
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Fig.  21 . Relative  reflectance  of  snow  and  clouds  as  a function  of 
wavelength  from  aircraft  observations.  (Figure  14  of  Valovcin 
11978).) 


can  be  distinguished  from  snow.  Figure  21  is  taken  from  the 
report  of  Valovcin  [1978].  The  reflectance  of  cirrus  clouds  is 
shown  to  be  greater  than  that  of  snow.  This  is  because  cirrus 
particles  have  a much  smaller  effective  mean  radius  than  do 
snow  particles  (perhaps  20  /xm  as  opposed  to  200  pm)  and 
the  albedo  is  higher  for  smaller  ice  grains  (cf.  Figure  1).  The 
fact  that  cumulus  clouds  have  even  higher  reflectance  is  due 
in  part  to  each  of  three  causes: 

1.  Cumulus  clouds  are  thicker  than  cimts. 

2.  Water  is  less  absorptive  than  ice  in  this  spectral 
region. 

3.  Cloud  water  droplets  are  smaller  than  cloud  ice  parti- 
cles. (The  cumulus  clouds  studied  were  all  water  droplet 
clouds,  even  at  their  tops.) 

These  aircraft  observations  have  led  to  the  employment  of 
a ‘snow  cloud  discriminator’  channel  at  X = 1.6  /xm  on  a 
Defense  Meteorological  Satellite  Program  satellite.  This 
channel  cannot  distinguish  between  thin  cloud  cover  over 
bare  ground  and  snow  cover,  but  these  can  be  distinguished 
in  the  visible  channel;  a bispectral  method  should  therefore 
be  quite  reliable  as  a snow  cloud  discriminator.  Tests  of  the 
utility  of  the  1.6-/um  channel  are  under  way  by  comparison 
with  results  from  human  analysis  of  satellite  imagery,  who 
use  pattern  recognition  methods  [Woronicz,  1981),  or  from 
ground-based  observations  of  snow  and  cloud  (R.  Barry, 
personal  communication,  1981). 

3.  Thermal  Infrared 

The  temperature  of  a snowpack  should  be  detectable  from 
the  thermal  infrared  emission,  provided  that  the  dependence 
of  emissivity  on  viewing  angle  (Figure  186)  and  the  atmo- 
spheric transmissivity  are  accounted  for. 

4.  Summary  of  Snow  Parameters  Detectable  in 
Solar,  Infrared,  and  Microwave  Spectra 

Table  2 offers  a rough  allocation  of  the  snow  parameters 
affecting  solar  spectral  albedo,  thermal  infrared  spectral 
emissivity,  and  microwave  spectral  emissivity.  The  micro- 
wave  emissivity  depends  on  three  additional  variables  which 
are  unimportant  for  the  shorter  waves;  (1)  liquid  water 
content,  because  the  for  ice  and  water  are  very  different 
in  the  microwave;  (2)  temperature,  because  the  mim  of  ice 
(and  of  water)  exhibits  a significant  temperature  dependence 
for  X & 0.2  mm;  and  (3)  density,  because  snow  grains  are  in 
each  others’  near  fields  at  these  wavelengths. 
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TABLE  2.  Parameters  Affecting  Albedo  and  Emissivity  of  Snow 


Visible 

Solar 

Albedo 

Near-IR 

Solar 

Albedo 

Thermal 

Infrared 

Emis- 

sivity 

Micro- 

wave 

Emis- 

sivity 

Grain  size 

(+) 

yes 

yes 

2enith  (or  nadir)  angle 

(+) 

yes 

yes 

yes 

Depth 

yes 

yes 

Contaminants 

yes 

Liquid  water  content 

yes 

Density 

yes 

Temperature 

yes 

present. 


No  attempt  is  made  here  to  subdivide  the  microwave 
region.  If  that  is  done,  one  may  find  that  the  emissivity  is 
insensitive  to  some  of  the  parameters  in  Table  1 at  certain 
wavelengths. 

The  thermal  infrared  emission  can  be  used  to  obtain  snow 
temperature.  But  a multiplicity  of  factors  are  seen  to  affect 
the  optical  properties  of  snow  in  other  spectral  regions.  In 
order  to  detect  individual  snow  parameters  unambiguously 
from  satellite,  one  must  therefore  examine  the  snow  at 
several  wavelengths  simultaneously. 

O.  Recommendations  for  Modeling 
and  Experimental  Work 

To  test  the  theoretical  radiative  transfer  models,  more 
visible  and  near-IR  measurements  are  needed  of  monochro- 
matic flux  extinction  and  of  the  zenith  angle  dependence  of 
monochromatic  albedo  in  situations  where  all  relevant  snow 
parameters  (see  Table  2)  are  measured  simultaneously.  The 
spectral  emissivity  in  the  thermal  IR  should  be  measured  as  a 
function  of  density  and  viewing  angle.  Reflectance  measure- 
ments for  unusual  snow  conditions  would  also  be  useful, 
such  as  snow  covered  with  fimspiegel  and  rapidly  melting 
snow  which  may  have  a high  liquid  water  content.  Special 
attention  should  be  devoted  to  characterizing  absorptive 
impurities  in  snow. 

For  remote  sensing  applications,  angularly  detailed  mea- 
surements of  the  bidirectional  reflectance  for  various  wave- 
lengths, grain  sizes,  and  surface  conditions  should  be  given 
high  priority. 

Topics  for  future  modeling  include  the  study  of  near-field 
effects,  nonsphericity,  and  surface  irregularity  and  the  pre- 
diction of  bidirectional  reflectance  rather  than  merely  fitting 
it  empirically. 

Note  added  in  proof.  Ackerman  and  Toon  [19811  have 
calculated  single  scattering  from  atmospheric  ammonium 
sulfate  particles  containing  soot  and  have  shown  that  the 
absorption  is  indeed  enhanced  by  putting  the  soot  particles 
inside  the  transparent  sulfate  particles.  To  reduce  the  aero- 
sol single-scattering  albedo  from  1.0  to  0.85  required  10-20% 
soot  (by  volume)  if  the  soot  was  present  as  separate 
particles  but  only  5%  if  the  soot  was  located  inside  the 
sulfate  particles.  Similar  results  would  be  obtained  for  soot 
particles  embedded  in  ice  grains,  as  was  speculated  in 
section  HI  above. 
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